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A A master integral 



1 Introduction 



In this work we present a systematic approach to the numerical evaluation of one-loop multi-leg 
Feynman diagrams. One may wonder why or what is new and in what respect there is something to be 
gained after the seminal work of 't Hooft and Veltman Q , where the analytic answer to the same problem 
was given. 

There are several reasons for attempting a systematization of the problem.^ First of all we want a fast 
and reliable way of computing one-loop multi-leg Feynman diagrams for all kinematical configurations, 
in particular around the (possibly) singular points of the diagrams which are determined by the solution 
of the corresponding Landau equations [^. For each graph we have a leading Landau singularity, usually 
called anomalous threshold for the triangle and higher functions, which for almost all practical situations 
lies outside the physical region. However, there are important examples where the external legs of the 
one- loop diagrams are not representing the in/out states of the process under examination and, for a 
general treatment, we must go beyond the class of 1 — > 2(3) or 2 — s- 2 processes, to a situation where 
classifying a priori the physical nature of the anomalous threshold is a tremendous task. 

Furthermore we have sub-leading singularities which are the leading ones for the reduced diagrams 
originating from the primary one Q . In principle, all of them are fully accounted for by the analytic answer 
but, in practice, the corresponding numerical evaluation of the special functions that form the answer can 
be rather unstable. For complicated processes, like full one-loop corrections to 2 — > 4(6) scatterings, we 
have a rich phase space at our disposal and any automatized calculation must be designed in such a way 
to foresee all possible sources of numerical problems without having to stop the full evaluation process. 

A particularly severe example is represented by the decomposition of pentagons into a sum of boxes p. 
For a real multi-scale problem these boxes are to be evaluated in regions where the analytic result requires 
cumbersome analytic continuations . Our goal is to derive a numerical answer without having to worry 
about kinematics. Even more important, our approach does not suffer from the introduction of complex 
masses which, instead, is a burdensome aspect in the analytic way. At most some care is needed when 
we take into account that complex masses are to be understood as complex poles lying on the second 
Riemann sheet [@]. 

Another well known drawback of the standard analytic approach Q is the appearance of negative 
powers of Gram's determinants in the reduction of tensor integrals, which represents a particular annoy- 
ance since their zeros do not correspond, under normal circumstances, to true singularities of Feynman 
diagrams. Clearly, in any approach and after a certain number of symbolic manipulations, certain de- 
nominators will show up in the answer. It seems therefore reasonable to design an approach where the 
zeros of these denominators correspond to true singularities of the diagrams. In any case, a new strategy 
for the so called problem of reduction of tensor integrals is highly desirable. 

Furthermore, in any realistic calculation infrared divergent (IR) virtual corrections will emerge, and 
we have found it very attractive to use the same tools that are needed to solve the problem of numerical 
evaluation of multi-scale diagrams for classifying and computing residues of infrared poles and IR finite 
parts of IR divergent diagrams. In this respect, the interplay between the special class of IR singularities 
and the general one of Landau singularities is essential. 

In this work we make use of a recent proposal for numerical evaluation of multi-loop integrals Q and 
extend the work of [|lO| ] to all one- loop diagrams (up to six external legs). The global strategy is based 
on the use of the Bernstein-Tkachov theorem which, away from the leading and sub-leading Landau 
singularities, allows us to cast any one-loop integrand in a form well suited for numerical treatment. As 

'^Alternative approaches, based on differential equation methods, have been originated in M. 
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it will be shown in Section ^ a unique factor Bq is associated with each one-loop diagram GJ^ For each 
family of one-loop diagrams we then consider three different cases: Bq ^ 0, Bq = but G regular, 
and Bq = with G singular. For the last case, the method is based on a particular combination of 
Mellin-Barnes |12| and sector decomposition ||l^ techniques, as a result of which we are able to write a 
Laurent expansion for G around the singularity. Alternatively, we have been able to derive new integral 
representations that are particularly suited for numerical treatment around the singularities. 

To a large extent, the use of complex masses (poles) represents a significant simplification of the whole 
approach. For physical (real) values of the Mandelstam invariants describing the process, none of the Bq 
factors can be zero, and a straightforward iteration of the BT procedure is all what we need. 

To return to our original question on why devoting such an effort to the evaluation of one-loop 
diagrams, we may add that the results of our work must be seen as one of the many building blocks 
towards the computation of physical observables at the two-loop level. Here our strategy has been already 
outlined in Q: full numerical analysis of the two-loop content. Clearly we will have to include the one- 
loop part, and it is rather obvious that the two pieces should be treated on equal footing. Furthermore, 
our main emphasis has been on a general treatment, therefore extending previous works on the same 
subject 0. 

We do not present numbers in this paper for a very simple reason: in all cases where the analytic 
answer is known, we have found perfect agreement with the one-loop library of TOPAZO The 
resulting comparison, basically a long list of identical numbers, is not very illuminating. 

The outline of the paper will be as follows: in Section ^ we introduce the definitions and discuss 
the Bernstein-Tkachov theorem, the Mellin-Barnes techniques, and provide a short discussion about the 
physical/unphysical role of the Landau singularities. In Section ^ the generalized two-point function 
is introduced, while in Section |^ we start the presentation of the one-loop, three-point, C-family, and 
in Section |5| we discuss its infrared divergent configurations. The four-point, D-family, is discussed in 
Section ^ and the corresponding infrared configurations are presented in Section |^ The five- and six-point 
functions are analyzed in Sections |M 1^ and 10. 



2 Prolegomena 

— The Bernstein-Tkachov theorem 

The Bernstein-Tkachov theorem [^] (hereafter BT) tells us that for any finite set of polynomials Vi{x), 
where x = (xi, . . . , Xk) is a vector of Feynman parameters, there exists an identity of the following form: 

V {x,d)llvr^^\x) = Bllvr^ix). (!) 

i i 

where P is a polynomial of x and dj = d/dxj] B and all coefficients of V are polynomials of /ij and of 
the coefficients of Vi{x). Any multi-loop Feynman diagram can be cast into the form of Eq.(|l|). Iterative 
applications of the BT functional relations, followed by integration by parts, allows us to express the 
integrand in parametric space as a combination of (polynomials)'^ x logarithms of the same polynomials, 
with any given integer w > 0, therefore achieving a result that is well suited for numerical integration. 
The B coefficients of Eq. (||) are connected to the leading Landau singularity of the corresponding diagram 
while repeated applications of Eq.(|I]), after integration by parts, will bring the sub-leading ones into the 
final answer. 

^This factor is proportional to the determinant of the coefficients in the system of the corresponding Landau equations; 
S. Uccirati, PhD dissertation (2002). 
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For general one-loop diagrams we have an explicit solution for the polynomial V which is due to 
F. V. Tkachov |11|. Any one-loop Feynman diagram G, irrespective of the number of vertices, can be 
expressed as 

G= f dxV-^{x), (2) 



where the integration region is Xj > 0, J2j Xj < 1, with j = 1, . . . , N — 1, and V{x) is a quadratic 
polynomial in x, 

V{x)=x^Hx + 2K^x + L. (3) 
The solution to the problem of determining the polynomial V is as follows: 

7? = l_i|__^ll^, X' = -K'H~\ B = L-K'H-'K, (4) 

where the matrix H is symmetric. In our conventions a multi-leg function is specified by the number of 
its internal lines, without regard to the actual number of external lines (see Fig. |2| for an illustration). 

- Mellin-Barnes transform 

Another tool that will be needed in the following is represented by the Mellin-Barnes splitting 



1 



+1 oo 



(Q + A)-° = — - / dsB{s,a-s)p^~'Q-', a > 0, (5) 
z vr z 



I CO 



where p = 1/A, B is the Euler beta- function and 

< Re s < a, | arg Q — arg A | < tt. (6) 

The i5 prescription is essential in deriving the correct analytic continuation. In all cases Q will be a 
quadric in n variables and A will be a constant, essentially a BT factor. Therefore Q is always understood 
as Q — i6 with (5 — > 0+ . Since | arg Q — arg A | < tt is required we will assign a small negative imaginary 
parts to both Q and A. 

- Anomalous threshold 

For a given one- loop diagram G, let us denote with Bq the corresponding BT factor of Eq.(2) which 
is a function of the internal masses and external momenta of G. At the leading Landau singularity of 
G, the so-called anomalous threshold (AT) we have Bq = 0. Conversely, Bq = is the condition 
to have a proper solution for the system of Landau equations corresponding to G (see Section |5.3| for 
additional discussions). Note that AT is not directly related through unitarity to physical processes (cut 
diagrams). Most of our work is devoted to analyze specific algorithms for the numerical evaluation of G 
around Bq = for arbitrary values of the internal masses and external momenta. However, for a large 
class of applications these parameters are bounded to the physical region TZ and the relevant question is 
whether or not AT £ TZ. A complete discussion of this issue is technically rather complicated and beyond 
the scope of our work; here we illustrate a simple case, namely the one-loop vertex of Fig. || with the 
following configuration, = m, 3 = —M"^ < and pi = —r. For real vectors p2,3 it follows that TZ 
is defined by either r > 4M^ (s-channel) or r < (f-channel). The diagram has an AT at 

r = r^y = 4M2 fl-^y iff M^>2m^. (7) 
V 4 / 



In our metric spacelike p implies positive p . Further p4 — ipo with po real for a physical four-momentum. 
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The kind of singularity depends on the value of M^: if < < 2 there is no AT, if < m2 < 4 m 
there is an unphysical AT (0 < r_4r < Arm?) and, finally, if > Arm? there is a physical AT at r^^ < 0. 

— Notations 

Finally, to keep our results as compact as possible, we introduce the following notations: 



dSn= / dxi / dx2 ••• / dxn, / dCn = / TT dXj, / dS2 = / dX2 

Jo Jo Jo J Jo J Jo Jo 



dxi. (8) 



Likewise we also need 



■l — Xl r-X-l+Xl—X-z j-Xn-l+Xn-l—Xn 

dSn{X) = I dxi I dx2 ••• / dxn, (9) 

X\ J—X2 J —Xji 



where X is an n-dimensional vector. 



3 Generalized two-point function 

One of the building blocks for the evaluation of all one-loop diagrams is the generalized two-point 
function, i.e. a self-energy graph with arbitrary powers for the two internal propagators. The generalized 




rm,2 

Figure 1: The two-point Green function, 
two-point function can be written in terms of the functions 

Bn{a) = [ dxx" \xB{x)-i6\~", (10) 

Jo L ^ 

where n is a non- negative integer and Xsix) is a quadratic polynomial whose coefficients depend on p^, 
the external squared four- momentum, and rm,i^2, the internal masses (see Fig. |l]). Let x± be the roots of 
the equation Xb — i ^ = 0- If A(— p^, , m^) = 0, where A(x, y, z) = x'^ + y'^ + z'^ — 2{xy + xz + yz) is 
the usual Kallen function, the integration contour may be pinched between the two singularities x± as 
(5 — > 0+. If this pinching occurs, the integral is singular with a branch point of the two-particle cut [^]. 
Consider 

Xb{x) = ax^ + bx + c = a{x - Xf + \, (11) 
where we have introduced the auxiliary quantities 

X = -^, A=-^, ^ = 1? -Aac = X{-p^,rmlml). (12) 
la Aa 

In the following we describe a solution for Bn{a) which allows to extract its behavior around A = 0. 
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We start our derivation by analyzing Bo{a). First we assume < X < 1 and obtain 

2 ^1 2 



So (a) 



V(-l)'Xi / dx{aXf + X-iSy^ = ^{-iy XiJi{a,X), 

1 = 1 



(13) 



where Xi = —X, X2 = 1 — X and where the i 5 prescription is equivalent to the replacement o, A — 
a,X — i6. In the following it must be understood that both a and A are given an infinitesimal imaginary 
part. Note that with this prescription we have 



a 

Next we split Ji into two pieces 



[-) -i sign (a) (9(--) 
a a ■ 



^ |l/2=;^l/2^-l/2_ 
a 



Jl{a 



poo poo 

A) = y dx{aXf x'^ + X)-" - J dx{aXfx^ + X)-'^, 



(14) 



(15) 



where the first one contains the divergent part for A — > 0, while the second is finite. In the latter we 
perform the transformation x x^^^^ and we obtain: 



Ji{a,X) 



poo -I pi 

/ dxiaXfx^ + Xy-- / (ixx°-3/2(aXf + Ax)-". 
Jo 2 Jo 



The first integral in Eq.(^) can be cast into the following form 

dx (aXf x2 + A)-° 



TTTT-. r -D - , a 



(16) 



(17) 



which is valid for | arg(A/a)| < vr and Re a > 1/2. Here B denotes the Euler beta-function. Collecting 
the two pieces we obtain 



X^/^-<^a-'/^B{^,a-^)-\Xi\ / dxx"~3/2 (^^f + Ax)- 



^o(«) - 2 2^ ^xT 

1=1 ' " 

= AV2-"a~i/2i3(-,a--)- - ^ \Xi\ C dxx^-'/^aXf + Xx)-", 

1=1 



(18) 



where we have used the relation \Xi\ = (—1)' X/, consequence of < X < 1. If instead X < or X > 1 
the first term vanishes and there is no divergence. In general we can write 



/ dxx"-^/^ {aXf + X : 
Jo 



\jT{aXf) "2i^i(a, a-^; a + ^; ^), 

a — 1/2 2 2 aXf 



(19) 



where 2F1 is the standard hyper geometric function [12|. Let us define zi = —X/{aXf); then, for a = 1 
we obtain 



1 3 

' 2 ' 2 ' 



1 



-1/2 1^1+ 



1/2 



while for a = n we will use 



1 1 ^ 

2Fi(n, n - - ; n + - ; zi) 



l-z 



1/2 ' 



(3/2)„_i d"-^ 1.3 
(l)„_i (l/2)„_i dz^i o. o'^^J- 



2 ' 2 



(20) 



(21) 
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where (a)„ = T (a + n) /T (a) is the Pochhammer symbol and F is the Euler gamma-function. Note that 
both integrals diverge for a — > 1/2 and that the divergent parts cancel in the sum. At a = 1/2 and with 
= aXf /X we have the following result: 



BoU = E (-1)' sign (Xi) a-V2 in [^^ + (1 + ^2)i/2 



1=1 



(22) 



In the same way we find that for the general case the result is 



dx x^' 



Xb{x) -iS 



-a 1 



2 n 



1=1 k=0 



^k/2+1/2-a ^-k/2-1/2 B(t±l,a - ^ 



i-iyx 

Jo 





\ 






1-1 




I dxx" 


Jo 





k+l 



(a Xf + \x) 



(23) 



The integral in the second term of Eq.(p^) can also be expressed through an hyper geometric function: 



Bn{a) 



dx x" 



Xb{x) -iS 



-a 1 



2 n 



iEE 

1=1 k=0 



I vn—k 



-lyx 



1l 

\Xi\ 



fc+i 



^k/2+1/2-a ^-k/2-1/2 B{'^,a - 



\Xl\ 



k+l 



{aXiY 



a-{k + l)/2 



k+l k-1 
2F1 [a, a — ; a — 



(24) 



Finally, a special case will be needed in dealing with multi-leg functions: let a = 1 + e/2 and consider 
the expansion around e = 0; for n = 0, we obtain 



i?o(l + e/2)=i?o(l) + 2^o(l) + 0(e') 



(25) 



1=1 



E 

n=0 



n + l/2 



?/'(n + 1) - V'(l) 



n + l/2 



ln{aXf-i6) (-A)", 



where ip is the Euler psi-function. More generally, for a non-negative integer value of n, we have 



Bnil + e/2)=Bn{l) + -BUl) + 0{e^) 



B: 



^ 2 n 



Xi 



\Xi\ 



k+l 



{b'l-k + &2;fe) , 



(26) 

(27) 
(28) 



-i)'x"-'= 

1 = 1 k=0 

where 6'^.^ and b^.j^ are the terms proportional to e/2 that are obtained by expanding around e = the 
first and second term in the square brackets in Eq.(p^), respectively: 



6;.^ = A^-/2-l/2^-fc/2-l/2B 



1-k l+k 



l-k 



^(1) -ln(A-i5) 
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V'(l)-ln(aXf 



k-l- 2m)2 
^5)]}. 



aXi 



{2 + {k-l - 2m) [V'(m + 1) 



(29) 



4 Three-point functions (C-family) 

We introduce scalar one-loop A^-point functions, with > 3, according to the following definition: 

If., 1 



iTT'^T{N -2) 



[q^ + mf] • • • [{q + Pi+ ■ ■ ■ + Pn-i)"^ + w 



(30) 



where d = 4 — e is the space-time dimensionality and Vm{x) is a quadratic polynomial in the Feynman 
parameters Xi, . . . , Xjv-i- In this section we analyze an arbitrary scalar three-point function Co = 5*3 (see 
Fig. ^), defined by 

Co= / dS2V~^-''''{xi,X2), (31) 



where y is a quadratic polynomial 



V{xi,X2) = axi + bx2 + cxiX2 + dxi + ex2 + f — id = x^Hx + 2K^x + L, 



(32) 



whose coefficients are related to the internal masses and the external momenta by the relations Hij 
— Pi ■ Pj, L = m\ and 



^'^ = \ (Pi ■ Pi + "i2 - ^i)^ K2 = ^iP ■ P -pi-pi+ml- ml), 



(33) 



with P = pi + p2. 





Figure 2: The one-loop, three-point Green function. The second diagram, although having 4 external lines, is 
included in the C-family of 3 internal lines. 



Let us define the usual BT factors^ as -B3 = L — K and X = — K. In the following we 

will present integral representations for Cq which are well suited for numerical integration and cover all 
options for the choice of external momenta. 

''Note, however, the different sign in the definition of X with respect to 
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4.1 Evaluation of Co for ^ 

If i?3 7^ we can rewrite Cq according to the standard BT procedure of Eq.(||). It is convenient to 
introduce special notations, Xq = 1, X3 = 0, and V{i i + 1) to denote contractions, i.e. 



y(oi) = y(i,xi), y(i2) = y(xi,xi), y(2 3) = y(xi,o). 

In this way we obtain a simple integral representation 

Co = + /" dxi[r dx2 In V{xi,X2) ^^+1) InFCz'TTi)] }. 



(34) 



(35) 



Note that a zero of the corresponding Gram's determinant is not a problem, actually it even decreases 
the number of terms. Indeed let G3 = detif and denote by the co-determinant of the element Hij. 
If G3 = we obtain 



1 ^ ~ 

Co(G3 = 0) = - — / dxJ2i^i- ^^Vii i + 1)> 

3 JO 



(36) 



where B-^ = b^/G^ and X = —AK, showing that, in this case, Co is a combination of three Bq func- 
tions [||]. 

4.2 Form factors in the C-family 

The most general C-family integral that we consider in any field theory is 



/n m 
dS2Q{xi,X2)V'^'^/'^{xi,X2), Q(xi,X2) = ^ ^ anmXiX^. 
n=l m=l 



The integral representation corresponding to Eq.(|35|) will be 

1 



C[Q(xi,X2)] = 

n=l m=l 

where the C-coefficients are 

1 



B-, 



^0 



(2 + n + m) (1 + m) ' 



Cl,n = (Xi -X2)x^+'" In V{xi,xi)+5mflX2X^ In y(xi,0) + {l-Xi)xT In y(l,xi) 



(37) 



(38) 



Cm,n = (m ^2 xi - (2 -I- n -I- m) xi X2 + nXi X2) In y(xi, X2). 
Examples of standard form factors of the C-family (see Q for the general definition) are 

BzCij = j dS2C^j lny(xi,X2) + J dSiY^ lnV{ff+l) + K^j, 



(39) 



(40) 
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with coefficients 





= ^(Xi-3a;i), 

= ^(^2-3x2), 


Cu = ^Xi2Xl, 

1 1 

C12 = 2 ^12^:1, 




1 

~ 2 
= 0, 


X2X1, c?i = -Xoi, iiTii 
1 

Cl2= 2^01^1' -f^l2 = - 


1 

"~3' 
1 

6' 


C2I 


= xi {2xi - Xi), 




! 


C21 = 


= -2^2a;i, C2i = --Xoi, 


K21 = 


^-22 


= Xi {2X2 - X2), 


1 ^ 2 V 

C22 = - 2 ^1 




'-22 - 


= 0' c^2 = -^^oia;?, -fC22 


1 

~ 12' 


^23 


= 2x1X2 - ^ (Xi X2 + X2X1), 










<^23 


= -^xlXi2, 


= 0, C°3 = -^ 


Xoi xi, 







wfiere we introduced Xij = Xi — Xj and Xq = 1. For C24 we directly obtain 

C24 = j^-^ / dS2 lny(xi,X2), ^ = --7-ln7r, (42) 

where 7 = 0.577216 • • • is the Euler constant and d = 4 — e is the space-time dimensionality. 

It is perhaps the right moment to summarize our strategy for tensor integrals. Clearly a decomposition 
in form factors is always possible and, within our approach, the form factors are not plagued by the 
appearance of inverse powers of Gram's determinants. However, this is not exactly the way a calculation 
should be organized since one has to profit from the strategy of the BT procedure by mapping all 
integrals, for all topologies, into just one integral. Therefore, there is no need to perform a pre-reduction 
in momentum space. Furthermore, in any realistic calculation, the denominators in Feynman diagrams 
contain scalar products, e.g. q ■ pi etc, and one has to benefit as much as possible from algebraic 
simplification. However, this is not always the case because each of the external lines, say in our original 
diagram, can split into other two (three) lines (we have in mind vector bosons that should be considered 
as unstable particles and, therefore, always be attached to fermionic currents). In this case momenta will 
appear that are not present in the propagators of the diagram and no simplification will occur. For this 
reason, and for the reader's convenience, we list all the form factors in the one-loop families. 



4.3 Co at O (e) 

A two-point function with the insertion of a counterterm, which is needed in a two-loop calculation, 
is formally identical to a one-loop three-point function with one zero external momentum. Its expression 
at O (e) is needed and this follows from the BT algorithm quite nicely, 

Co{d) = Co (4) + e C^') + O (e^) , (43) 

Ci'^ = -\ J dS2 [lnF(xi, X2) + ^ ln2 F(xi, X2) 
1 r 2 

i=0 

showing that it contains at most squares of logarithms of quadratic forms. 



10 



4.4 Evaluation of Co when Bs ^ 



When we are around the region where B3 = 0, Eq.(pj5|) cannot be used and an ahernative derivation 
is needed. Actuahy there are two sub-cases to be discussed starting from the observation that Co can be 
written as 



Co 



dSo 



(x - X f H{x-X)+B3 



(45) 



If S3 = the integral of Eq.(H5|) is singular or regular depending on the values of Xi^2- 



4.4.1 The regular case: method I 
If S3 = 0, but the condition 



< X2 < Xi < 1 



(46) 



is not fulfilled, there is no real singularity inside the integration domain X3; still we cannot apply the 
standard BT techniques of Eq. (||) . The integral in Eq. (^) is well defined through distortion of the contour 
integration and there is no pinch inside X3. We use this regularity to write the following BT relations, 
corresponding to S3 7^ and S3 = 0: 



1 



2(/i + l) 
1 



{X-X)d^ l/^+l(xi,X2) = S3l^'^(xi,X2), 

-i+l, 



2(/i + 

with V = Vq + B3. Subtracting the two equations, we obtain: 

1 . .■\Vt'+Hxi,X2)-VQ^'-\xi,X2) 



(47) 



V^{xi,X2) 



1 



2{fi + r 



(x - X) 



S3 



(48) 



Note that this relation holds even in the case S3 ~ 0. Since the point of coordinates Xj = Xi is outside 
X3, we can perform an integration by parts on both V and Vq. The result is the following: 



dS2ln^-l I dS, Y: (^^ - In ^^'ill'^ 



i=0 



(49) 



All tensor integrals can be treated according to the same procedure. Note that the logarithms appearing 
in Eq.(^9|) are of the form ln(l + x) with small |x| and the corresponding values must retain full relative 
precision, e.g. they must be computed with routines based on Chebyshev expansion. 



4.4.2 The regular case: method II 



Even if the method of Section 4.4.1 gives a compact result for the regular case, we derive another 
expression for numerical checking of our results. Since Co is not singular in the limit S3 — > we can 
use a Taylor expansion of the integrand of Eq.(|3l|) followed by a BT-algorithm corresponding to S3 = 0. 
Therefore we have: 



00 „ 00 

Co = E(-^3)" / dS2{X){x'Hx)-^-'~^/^ = Y,C{n + l){-B,r, 

n=0 n=0 



(50) 
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where we kept e 7^ only where strictly needed and where we used Eq.(^). By using: 

1 



dS2{X) 



1 + 



2n + e 



X dr 



(51) 



We observe that integration-by-parts is allowed since = is outside the integration domain and derive, 
after changing x into either x — Xi 01 x — X2, 



C{n)= / dS2{X)V"'-'/\x^,X2) 



2n 



1 r 2 

/ dSiY.{X,-X,+^)V-^~^/\Tii + l), (52) 

i=0 



where we employed Xq = 1 and X-^ = 0, and we introduced shifted arguments, i.e. the three relevant 
quadratic forms in x are 

V{t5\) = Vc = V{1- Xi,x - X2), V{T r2) = Va = V{x - Xi,x - X2), 

V{T22,) = Vb = V{x-Xi,-X2), (53) 

for which we write Vi = hix"^ + 2kiX + li and introduce the standard BT factor, bi,i = a,b, c, and also the 
BT co-factor, Xi,i = a,b, c. As long as all the hi are different from zero, the corresponding BT relations 
hold. For the first term in the expansion we obtain: 



^(1) = ^ Yl it ^0;^(l)+/ dxCi4l,x) 



XI 



i=a,b,c 



C0;i(l) 



X, ln2 ^^(0) + (1 - Xi) In^ V^il) , Ci;,(l,2;) = InV'(x) + - In^ Viix), 



where the leading coefficients are: 

Xa = Xi - X2, Xl; = X2, 

whereas for n > 1 we get 



X^ — 1 - Xi, 



C{n) 



4(n - 1) 



1 _ ^ 
1)2 2^ 



bi 



Co-iin) + / dxCi-i{n,x) 



Co-i{n) = (1 - X,) Vr^+'{1) + X, Vr^-^'{0), Ci,i{l, x) = (2n - 3) V, 
Now we can sum the series according to the weU-known results: 

CxD „ 00 „ 



-n+l. 



X 



n=l 



n=l 



where Li2 (x) is the standard dilogarithm |2C], and obtain a novel representation for Cq: 



i=a,o,c 



Co-i + dx Ci-. 



(54) 
(55) 

(56) 

(57) 
(58) 

(59) 
(60) 



Co-', 



2 Lis 



_B3_ 

^^(0) 



+ ln2y,(0) 



C,., = \n[Vi{x) + B,]+-U2 



4(1-^0 



2U2 



(61) 
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The above derivation of the coefficients in the Taylor expansion, Eq.(l52|), assumes that ah sub-leading 
BT factors are non zero. C{n) of Eq.(52) is a combination of generalized two-point functions of Section |3|; 
if one of the 6j ~ we still have to distinguish between two possibilities: 

a) bi = but -hi < h < {hi > 0) is not fulfilled. 

b) bi = and -hi <ki<0{hi> 0) is fulfilled. 

In both cases we can apply the results of Section 0. So let us denote each term of Eq.(^) with 



1 



dxgiV, " '^^(x) 



2n-2+e 

where gi = —Xi + Xj+i. The corresponding contribution to Co will be 



(62) 



(63) 



ra=0 



In evaluating the sum it is more convenient to isolate the leading term CQ(n = 0) and consider the rest 
separately. The latter becomes 

{-l)'Xl ^ (-B3 



n=l ^ 1=1 ' * ' n=l 



\h-"^"'K'"B{\,n+\)-\Xl\ I dxx'^ih,xf + b.x 



n 



2\-n-l 



(64) 



where Xu = —Xi and X2i = 1 — Xi. The second series can be easily summed up, while for the first one 
we have: 



E 

n=l 



-^3 A)" ^ , 1^ _„i/2 f - (n+l/2)r(n+l/2) 

{n + lYT{n+l) 



n 



n=0 



/■I 1 

= — TTKj / dx {1 -\ — In x) {I + KiX) 
Jo 2 

where we introduced Ki = B^/bi. Finally we obtain: 



\n+l 



-1/2 



(65) 



J rest 



-lyxu 



t \Xh 



dx 



TT Ki X + x\nx 



{bihiYl^ (1+K,X2)V2 



\Xu\{hiXt^ + bix')-' In 



2 1 h^Xl + {bi + B^)x'^ 



hi Xf. + bi x2 



(66) 



As we argued for the generalized two-point function, the divergent term for bi ^ cancels out if the 
corresponding Xi is outside the interval [0, 1]. For the leading term, after an e expansion, we have: 



C^(n = 0)= l\xg,Vr\x) [i - i 



1 ^ 

\nVA, V / dxgiVr^ = Q. 



(67) 
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Note that the divergent term does not contribute to Cq. In order to calculate the term Cg(n = 0), we 
need the relevant part of C(l), i.e. the one corresponding to those values of i for which bi ~ 0. The latter 
contribution can be obtained by observing that, symbolically, 

eC{l) = Y,Bo{l + e/2y, (68) 

using Eq.(|2^) we expand the relevant generalized two-point function and derive an expression for C(l) by 
equating the terms proportional to e. 

Also for the form factors we can derive a similar expansion. Consider the auxiliary function 



^,-n-l-e/2 



n=0 i=l 
oo ni 712 / \ / \ 

E E E (?) (?) xr-^^xr-'^c{n+iAM){-B,T, 

^—n t.—n /„_n \ 1 / \ 2 / 



n=0 /i=0 /2=0 

and the expansion is fully specified by 

1 



C{n,ni,n2) 



dx { [{x - Xir - X^T (^1 - ^2) K"" 



2 — 2 n + ni + n2 

- (x - (-X2)"2+^ V^'' + (x - X2T^ (1 - y-"}. (69) 

Combinations of these functions account for all form factors. 



4.4.3 The singular case 

If condition Eq.(|4^) is fulfilled then there is a pinch inside X3 and the integral in Eq.(^) is singular. In 
order to minimize the number of terms in the solution we write 



Co= dS2V-^-'/^{xi,X2) = / dC2V-^''/\xi,X2) - / dxi / dxs F"^"'/^(xi, xs) 

J J Jo Jxi 

= J dC2V-'-'/\xuX2)- J dS2y-i-^/'(x2,xi) = Cr^'-^-C7o^°"'P. (70) 

Since the point Xi = Xi is now internal to the integration domain the complementary Cq function will 
be regular and can be computed according to the results of Section 4.4. 1| and 4.4.2 . A special case is 
represented by < X2 = Xi < 1 where both the original Cq and the complementary one are singular. 
However, due to the symmetry of H, we easily obtain 

Cq = C™ = ^ C™^ for 0<X2 = Xi<l. (71) 

^square ^j^^ integral over [0, 1]^ which we rewrite by changing variables, x'^ = Xi — Xi, as 



^square 



1=1 

where 



^a,A J dC2[Q^{xl,X2)+B3] ' (72) 

i=l 

ai = "2 = 1 - Xi, 03 = a4 = Xi, /3i = /33 = 1 - X2, (32 = Pa = X2, (73) 
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and where the new quadrics are defined by 

Qi = Q((l - Xi)xi , (1 - X2)x2), Q2 = Q((l - , -X2X2), 

Qs = Q{-XiXi , (1 - X2)X2), Qi = Q{-XiXi , -X2X2), 

with Q = x^ H X. In general we define 

Qi{xi, X2) = Aix\ + Bi xl + Ci xi X2. 



(74) 



(75) 



In order the derive a Laurent expansion around S3 = we introduce = and perform a Mehin- 

Barnes splitting as described in Section 0, 



^square 
^0 



Y^iPr dsB{s,l-s)p\~'Ci{s), Ciis)= dC2Qi%Xl,X2). 

Let us consider in detail the Cj-functions. We use a simple sector decomposition [13| to obtain 



(76) 



dS2+ / dS2 



Qi'{xi,X2)= / dC2xl {Ai + CiX2 + Bixi 



r 1 ^ 

+ j dC2xl-^'iB, + CiXi + Aixl)-' = jj:^—^J2Ci,,{s). 



(77) 



For each of the Cj^-functions we have a reduced quadratic form in one variable. Let us postpone for a 
moment the problem of their evaluation. From Eg. ([7^) we obtain 



i=l j=l J -too 

— I dx {flij X 2 k'lj X -\- lij^ 5 
^0 

and < Res < 1. Suppose that the corresponding BT factors 



r(s)r(i^ i_ 

P3 ^ijisj, 



lij - kfJhij = af Pf det{H)/h 



(78) 



(79) 



are not zero and that we are interested in the region of large |p3|. Then we close the integration contour 
over the right-hand complex half-plane at infinity. The poles are at s = 1 (double) and at s = k + 1 (single) 
where k > 1 is an integer. For Cij we use 



Cij(s) 



dx 



'J'ij Jo 



1 + 



1 x — Xii d 



2 s-1 dx 



Qli'ix), 



where Xij is the BT co- factor for Qij. In the hmit s — > 1 we obtain 
C,,{s) = C,,(l) + {s- l)C^(l) + 0{{s- If) , 



(80) 



Oij Jo 
Ojj Jo 



dx 



1 - - (x - Xij) — In Qij{x) 
' dx 



dx 



In Qij{x) + - (x - Xij) ^ In^ Qij{x) 



(81) 
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Therefore the residue of the double pole at s = 1 is 

Tlij \s=i= Cij{l) lnp3 - C'ij{l). (82) 
For the single poles at,s = fc + l,A;>lwe find that the residues are 

\s=k+i= -^-^J- Bo{k + 1) ^ (83) 



where BQ(k + l) is a generalized two-point function discussed in Section]^. With the help of Eqs.(^)-(p3|) 
we obtain the full Laurent expansion. 

Also for the form factors we can use the same strategy and write the auxiliary functions C(ni, 712) as 
a combinations of Mellin-Barnes integrals 



+^°° , r(s) r(i-s) , , 

1 - s + (ni + n2)/2 ^3 ^JK J, ), 



Cij{nj, s) = dx x^^ {hij + 2 kij x + li 
Jo 



(84) 



The relevant poles are ats = A; + l, fc>0 and there are two possibilities: if ni + 722 = 2n then there 
is a double pole for s = n + 1, otherwise, if ni + 71-2 = 2n + 1, there is one additional, simple, pole at 
s = n + 3/2. 

However, the secondary BT factors, bij, could be near to zero. Note that the bij coefficients are 
sub-leading of a second type since bij = should not be confused, in general, with the condition for 
the occurence of a Landau singularity associated with a reduced diagram. Indeed, following Eq.([7^) and 
Eq.(l79|), they correspond to having = X2 < Xi < 1 or < X2 < Xi = 1 in the pinch singularity 
and it is enough to select rrii = m, Mi in Eq.(|33[) to see that there is no correspondence with any of the 
two-particle cuts. From Eq.(^) and Eq.(^) it is easily seen that for G3 = dei{H) ~ we hace bij ~ G3 ^ 
so that G3 = does not imply bij = 0. 

bij = alone is not yet a sign of instability because Cij will show a pinch singularity in [0, 1] only if 
hij > and — hij < kij < 0. If we are in this region for some pair ij then it is better to use another 
expansion algorithm. From Eq.(^) we see that the total result is a sum of terms with the following 
structure: 

1 f'"'^ .„r(.)r(i-.) 



pF{p) = —- / ds ' ' \ '-f{s)p 

2 IT I . 1.:^. S - 1 



where a = k/h. With = 1/b we write 



1 00 

l+a 



dx{hx' + b)-\ (85) 



/(^) = pIl {fj'^ - £} dx (1 + PsLhx')-'. (86) 



Introducing a second Mellin-Barnes splitting gives 

2 01 o-l-joO 



[s - 1) (t - 1/2) \ p J 
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with Pi = a and /32 = 1 + This representation is vahd in the vertical strip < Re t < Re s < 1 and 
Ret < 1/2. Since we are interested in the region where both the leading B and the sub-leading b are 
approaching zero, we close the t-contour over the right-hand complex half-plane at infinity with poles at 
t = 1/2 and t = s + k,k > 0. We choose Res > 1/2 and obtain 



pPn 1 



n=l 



2 2iri 



+ 1 oo 



ds 



vr 



1/2 p 



I CO 

oo 



r(i-g) /pslV 

s-l \ p ) 



k=0 



k\ s + k-1/2 



T is naturally split into two parts: 

2 



-1) 



Hp. ps.) = E (-1)" l-^i + E ^ -^2 

n=l A:=0 
2 



(n,fc)' 



E (-1)" ^ [-^/^ ih ai ps.r'i^ Ti,n + E ^ Hn, 



(89) 



n=l fc=0 

For the first part in Eq.([89|) we introduce the new variable k = PslI P and derive the following coefficient: 

1 



27r? 



+ioo _ s) _ 1/2) 

as : K" . 

—i oo 



S - 1 



(90) 



We close the contour over the left-hand complex half-plane at infinity, with poles ats = l/2 — A:,A;>0, 
obtaining 



l,n 



g (-1)^ r (fc + 1/2) 



.1/2 



E 



r2(A; + l/2) {-k) 



k\ k + l/2 ^ r(/!:-F3/2) k\ 

fc=0 ' fc=o ^ ' / 



I/O ^ ,1 1 3 1, „ I/O , I-F(I-I-k) 
= -2 vr 2Fi(- ,-;-;--) = -2 vr^^ ^ in 

As a consequence we have 



1/2 



in) 



■ (R\ /^SL^l/2i P^'^ + (P + Psl)^/^ 

■ Sign {(3n) vr (— ) ^ In -— . 

pL 



(91) 



(92) 



For the second part in Eq.(^) the poles are at s = / -|- 1, / > where / = gives a double pole. We obtain 



•^2,n,, 



r(fc + i) 

k + l/2 



{hdir'-'p^'psi 



1 



k + l/2 



+ -1 + il,{k + l)-\n{hplp) 



_f{-Vl T{k + l + l) 
^ II l{k + l + l/2)^'^'' 



P PSL ■ 



(93) 



Substituting back into Eq. (|76|) gives the requested result. 

There is another situation that requires particular care: consider a quadratic form hx'^ + 2k x + l and 
suppose that b ^ but that / or that, equivalently h + 2k + I ^ 0. The solution that amounts to 
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iterate the BT procedure for f{k) is unstable because of the surface terms in the integration by parts. In 
this case we take the Mellin-Barnes transform of the integrand and write 



■ — : / dtB{t, s-t)L'~^ / dxx~\2k + hx 

'^^ J-ioo Jo 



1 



27ri 



too jO 

dtB{t, s-t)L'-^ ^ ' 2Fiit, 2-t; 

-7 nri 1 ^ 



2k' 



(94) 



where L = 1/1 and < Ret < Res < 1. Therefore we obtain 

P 



Hp,L)^pF{p) 



(27ri)2 

{2kLr\Fi{t, l-t-2-t- 



h 

2k 



(95) 



For large | L | we close the t— integration contour over the right-hand complex half-plane at infinity with 
poles at t = 1 and t = s + n. We get 



Hp,l) 



p 



(27ri)2 

X {2k L)-'- 

For the first term we write 
1 



■ioc S-l \pJ 

2-^1 (s n , 1 - s - n ; 2 



°° ('_-|\ra 

r(s-l) {2kL)-^ + Y^ ^^T{s + n-l) 



n=0 



n \ 



2k' 



27ri 



+*°° ^^^ r(i-g) r(s-i) /L 

■ioo s-l \p 



n=l 



P- 



(96) 



(97) 



where we closed the s— integration contour over the left-hand complex half-plane at infinity, with poles 
at s = 1 — n,n > 1. For the second term we need special cases of the hypergeometric function [12|: 

^2n+l- 



2-Fi(-n , l+n; 2 + n; z) 



1 



(2 + n) 



1 - 2\"+i d" 



dz" L 1 - z 



(98) 



4.5 A new integral representation for Co 



It is very important to derive our results, for all kinematical regions, in more than one form in order 
to be able to provide reliable cross-checks. For this reason we derived another representation for Co, to 
be used when the BT coefficient -63 vanishes and the integral is singular. Starting from Eq.([f5[) a sector 
decomposition of the square [0, 1]^ gives 

4 



Cr"'"' = ^a,A dS2+ dS2 



i=l 
4 2 



Qi{xi,X2) + B3 

l-e/2 



-l-e/2 



(99) 



= ^^aift / dC2X2 xlQij{xi) + Bs 
i=i j=i 

where we introduced Qii{xi) = Qi{xi,l) and Qi2ixi) = Qi{l,xi). They have the following form: 
Qiji^i) = hij xf + 2kijXi + lij. Note that, according to the i6 prescription, the following replace- 
ment must always be understood: QijjB^ — > Qij,B3 — id. Now we use the results of Appendix A: the 
coefficients in the corresponding master integral are: 

a = Qij{xi), b = B3, a = l + e/2, (3 = 1. (100) 
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dxi I ^ ^3 Q 



Therefore we obtain: 

i=i j=i 

Qij{xi,X2) = Qij{xi) + Bsxl- 
different from zero we apply the following BT relations: 



dx2X^^"^'Q-j^ '^'^{xi,X2)\, 



(101) 



If the factors bij = kj — kfj/hij = a^Pfdet{H)/hij, corresponding to the quadratic forms Qij{xi), are 



1 



1 



h- ■ 



dx hiQij{x) 



Q.-.(.,.,) = _|l-_ 

Let us introduce new auxiliary variables: Qji 
result: 



C square 




xi + ^]dx^+ X2 



Qtj{xi,X2) 



(102) 



(103) 



kij/hij and Qj2 = 1 + kij/hij to obtain a compact 

(104) 



4 2 

E\ ^ Olj Pi ^ 



Cij = 2\nB^ 



dxi \nQij{xi) + Cij2 InQij(l) + Ciji InQij(O) - 2 



1 

dxi I dxq 







lnQij{xi,X2] 



dxo 



Cm 



X2 

In Qij{l,X2) 

X2 



+ / dxi 



\v? Qij{xi) + A\nQij{xi,l) 



ijl 



In Qij{Q,X2] 

X2 



\n^Qij{l)-Qn ln2Q*j(0). 

(105) 



All the Qij and -B3 are given a —i5. The '+ '-distribution is defined, as usual, by its action on a generic 
test function f{x): 

7(^) 



(106) 



If instead some hij approaches zero we must repeat the procedure. This can be achieved if we start from 
Eq. (llOl| ) and change variables according to xi ^ xi — kij/hij. Then we decompose the xi integration 
interval and remap each piece into the unit square; the result is 



4 2 2 



^square 



1 „_ 



^ ^ ^ aj /3i Ciji / dxil- ^'"^ hij (fji x\ + hj - ib 
i=\ j=\ i=\ \^ 

dX2 X2^^^ 



f 

Jo 



hij Ciji ^1 +B3XI + bij 



-l-e/2 



(107) 



where the replacement bij, hij — > bij, hij — 16 is also understood. Now we apply the result of Appendix A 
to both integrals, deriving 



i=i j=i 1=1 L 



— (h--c\)-^/^ 5-^/2 5-1/2 

2g V"ij ^iji) -D3 "ij 
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1 „ /I 1 e\ 



-Br'" 

1 



2' 2 2 

1 

dxi 







/ dX2 X2 




Jo 





-^l/2-e/2 







+ / dxi I dx2 xl x^ hij + ( i?3 X2 + hj ) Xi 



-l-e/2 



(108) 



The first three terms in Eq.( |108[ ) are divergent for S3 ^ bij 0. The second integral depends only on 
the ratio 

^ (109) 



bij — i6 ' 



so that, when we apply the master formula of Appendix A, the result is 



dX2 X2 ^^'^ 



B3 xi + bi 



-l/2-e/2 1 



-1/2 „-e/2 



-b'r Br' B 



2 



e 1 



2' 2 



dX2 



bij X2 + -B3 



-l/2-e/2 



• (110) 



Note that the first term in Eq.(llO) cancels the first term of Eq.( |108| ): 



i=i j=i 1=1 

1 „ /I 1 e 

'2 



C^7 + bij x\ 



-1 



+ 2'2 + 2) ^S'^"'^' 



dX2 



bij X2 + -B3 



+ / dxi I dx2 x\ X2 /ijj + (-63X2 + 6ij ) x{ 



-l/2-e/2 
-l-e/2 



(111) 



Using this intermediate result we can apply the BT procedure to increment the power of the quadratic 
forms in the first and in the last integrals while the second one can be calculated analytically. Finally a 
Laurent expansion around e = is performed giving: 

4 2 2 



-square 




1 



i=i j=i 1=1 



^square 



+ 



'•13 sijl 

1 

h../-2 
"-IJ '^ijl 



1 + (1 + K,,)^/' 



.1/2 



lnS3 1-- lnPi,,(l,0) +- ln^Pi,7(l,0) 



^square 

C square 
1 



1 



'hi I-,., 



- In Bs In Piji{x, 0) - -In ' 



1 /lnPij7(xi,a;2) 



a:;2 



where we introduced the new quadratic form: 

Piji{xi,X2) = hij Cfji + bij xl + Bsxl- id. 
Note that the poles in e cancel out in Eg . (11121) , as expected. 



(113) 



(114) 
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5 Infrared divergent Co 



Before discussing the general classification of IR divergent configurations for the three-point function, 
we consider a specific example, namely pf = —m\,p\ = —m^ and m2 = (see Fig. §). We obtain 

C^"" = Co(Q', -ml -ml ; mi, 0, mg) = J dC2 x^^'^ X'^~''\x2) = --^ j\x x'^'''\x), (115) 
where we introduced 

X{x) = -Q"^ + {Q"^ + ml - ml) X + ml, Q^ = {p^+p2)^. (116) 
In this case the sector decomposition is trivial. The BT factor and co-factor associated to the quadratic 

1 



form of Eq.(116) are 



5, 



u ^2 2 2^ Y ^ Q^ + ml-ml 



(117) 



and -Bth = at the normal(pseudo)-threshold of Cq. In the following we discuss the various cases under 



which Eq.(115) can be classified. 



5.1 Bth ^ 

As long as B^Yi / it is possible to apply the standard BT procedure; we derive 



1 



5th e 



dx 



1 + ^ (x - X) ^] [l - ^ Inx(x) + J ln2 x(x) + O (e^) 



where the infrared residue and the finite part are 



(118) 



Cr = -l + l{l-X) lnx(l) + ^Xlnx(0)-^ dx\nx{x), 

'^0 



(1 - X) ln2 x(l) + lxin^ x(0) + \ ^x 



Inx(x) + - ln^x(a;) 



(119) 



5.2 Infrared configuration and Sth ~ 

If -Bth ~ then we introduce h = —Q^ and write 



'^0 



- / dx{hx^ +B,^)-^-'/^. 
e J-x 



(120) 



and consider three sub-cases: from Eq.( |117D we see that -Bth = only if < 0, i.e. for —Q^ = (miztms)^. 
Note that, in this case, we have X = m\l{m\ ± ms). 
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5.2.1 0<X <1 



Suppose that < X < 1 so that the configuration is actuahy singular at -Bth = 0. This corresponds to the 
normal threshold = —{mi + m^)'^. Introducing a Mellin-Barnes splitting, as described in Section |2|, 
and performing the x-integration gives 



1 1 



2e 27r« 
X {h — id) 



+ ioo r(-s) r(s + l + e/2) 

as 



r(l + e/2) (s + l/2 + e/2) 



(121) 



with —1 — e/2 < Res < —1/2 — e/2. Note that we explicitly indicated the infinitesimal imaginary parts. 
Since we are interested in the behavior for Bth — > the s contour can be closed over the right-hand 
complex half-plane at infinity, with poles at s = —1/2 — e/2 and at s = k,k > 0. The residue at 
s = -1/2 - e/2 is 



At s = k we have instead 



R{s = k) = j-^h-''-'B^,^ E ^(«)«''" 



a=XA-X 



R{a) = { 



2k 



k + l/ 

1 1 



-k + l/ 

Finally, we present the last two sub-cases. 



^ +ln(/i-i5) +2 lna + V(l) -V'l^ + l)!)- 
VIZ J J 



(122) 



(123) 



5.2.2 X <{) or X>1 

In these cases, corresponding to the pseudo-threshold = —{mi — m^)"^, we will write 



r<iR 



1 r 

e 



9{-x) - e{x - 1] 



dx — dx 
Jo 



{hx^ + B,i,)-^-^/\ 



Suppose that X < 0, after the Mellin-Barnes splitting we obtain 



1 1 



2e 27ri 



+ ^oo Y{-s)T{s + l + e/2) 

CIS 

r (1 + e/2) (s + 1/2 + e/2) 



(/,_i5)— 1-/2 {Bii,-i6r \{l-X)-^ 



-2s-e 



-X) 



There is no pole at s = —1/2 — e/2 but only poles at s = k,k > with residues 



R{s = k) 



(-1)' /^-fc-lRfc 



k + l/2 



th 



{l-X)-^-'^^R{l-X)- \X 



-l-2k 



' R{\x\), 



with a similar result if X > 1. i? is defined in Eg. ( 123 ). 



(124) 



(125) 



(126) 
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5.3 IR configurations and Landau equations 

The study of the Landau equations for a given diagram is the most elegant way to classify its infrared 
divergent configurations. Again we refer to the one-loop vertex for illustration. The Landau equations 
are 

ai {q"^ + ml) = 0, a2 {{q + Pi)^ + m^) = 0, {{q + pi + p2)^ + mj) = 0, (127) 
ai q,i + a2{q + pi)^i + as{q + pi+ P2)n = 0. (128) 

A solution will be the leading singularity if it corresponds to ai ^ 0,V«. Next we multiply Eq.( |128D by 
q,pi and p2 respectively to obtain an homogeneous system of three equations where, moreover, we use 
q^ = —m\ etc, from Eq.( |127| ). A necessary and sufficient condition to have a proper solution, i.e. not all 
the ai = 0, requires the determinant of coefficients to be zero, thus fixing a relation between internal and 
external masses. Any configuration that satisfy this constraint is a Landau singularity for the diagram 
which, however, does not necessarily imply that the diagram itself diverges at that configuration. 

Suppose that we consider the following case: Pi2~ ~"^^5"^i,3 = ^^^d m2 = 0. It is easily found 
that this configuration is a Landau singularity. The question is however of which kind. Let us insert the 
above values into the homogeneous system, what we obtain is 

m^ai + {]^Q'^ + m^)as = Q, (^ + m^) ai + ag = 0. (129) 

First of all we observe that it is not either ai = or = —mf, etc; it can be both. Secondly, our 
configuration, where = (pi + ^2)^ is unconstrained, is a singularity. Finally there is a special case 
of the general configuration discussed (Q^ free) which is even more singular giving, in the infrared case, 
the true leading Landau singularity. To have 01^3 7^ we must require = — 4m^ which gives, in 
the annihilation channel, the well-known threshold singularity on top of the infrared one. This condition 
emerges also from the following argument: inserting 2 = —m'^ in the condition to have a proper solution 
one obtains 

mlQ'^ = 0, and (m^ - - 4m^) = 0, (130) 

In a certain sense the constraint = — 4 m^ is buried inside the anomalous threshold condition. In other 
words, that all the propagators in a diagram are on-shell and that the consistency relation is satisfied 
does not necessarily imply that all ai are different from zero: the infrared case is a clear example. Note 
the presence of a potential singularity also at = which, however, is not physical, i.e. not on the 
physical Riemann sheet. The latter fact can be seen by inspecting the explicit analytic result, 

Our strategy in the general classification of infrared divergent configurations, diagram by diagram, will be 
to assume a certain number of zero internal masses with at least one unconstrained external momentum. 
Then we fix the remaining parameters to satisfy the consistency relation for the Landau equations. 
Finally, we return to the original set of Landau equations and look for additional constraints that are 
necessary in order to arrive to the true leading singularity. The presence of a threshold-like singularity 
on top of the infrared poles is the sign that, after extracting these poles, we still have complications for 
the residual integrations that cannot be solved with naive methods. 

For Co and arbitrary masses we look for a solution of the homogeneous system of the three Landau 
equations in a situation where is not constrained (due to the symmetry of the problem it is sufficient 
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to consider Q^). The condition to have a proper solution is a quadratic equation in with m| as the 
coefficient of therefore we require = 0, after which we obtain the fohowing condition: 



{p'i + mj) {pI + ml) - {p'i + mi) - {p2 + mg) m{ {pi + m'^) - m'i {pf + m'() 



,2 / 2 



,2 / 2 



0, 



(132) 



requiring pi = —mf and = ~"^3) i-e- the configuration that we akeady discussed. If we substitute 
back into the Landau equations we obtain 



2 ai + (Q2 +ml + m|) = 0, 
(Q^ + mf + mj) ai + 2 mi Q3 = 0, 

so that the leading singularity, i.e. all a's different from zero, is obtained for 

Xi-Q\mlml)=0 = 8^ = 0, 

which clarifies once more the physical meaning of the BT factor. 



(133) 
(134) 



6 Four-point functions (D-family) 



Pi 



mi 



7712 



P2 



PA 



7774 



7773 



P3 



Figure 3: The one-loop, four-point Green function of Eq.( |l35| ). All momenta are flowing inwards 
The generic, scalar, box (see Fig. |3|) is given by 

Dq = j dS:iV-'^-^/'^{xi,X2,x:i), V{x) = H x + 2 x + L, 
where Hij = — pi ■ pj, L = m\ and where 



Ki = {ml -ml-pi -pi), K2 = (7772 - ml- 2pi ■ P2 - P2 ■P2), 



K3 



1 



[ml - ml- 2pi ■ p3 - 2p2 ■ P3 - Ps ■ Pz)- 



We are now in a position to write the answer for the case -B4 ^ 0. 

6.1 Evaluation of Dq when B4 7^ 

If we are away from B4 = L — H^^ K = we can apply a BT iteration obtaining 

3 

BiDo = -^ J dS3V-^-'/\xi,X2,X3) + ^ J dS2^{X,-Xi+i)V-'~'/\frri) 

i=0 



(135) 



(136) 



(137) 
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where X = —KH and, moreover, Xq = 1,^4 = 0. Contractions are defined as before with 



1 = (l,xi,X2), 3 4 = (xi, j;2,0). 



(138) 



etc. Therefore DQ[d = 4) is a combination of a non-contracted term which can be seen as DQ{d = 6) 
and of four pinches which can be represented as the triangles arising when we shrink a hne of the box to 
a point. The latter are three-point functions that can be treated according to the results of Section ^. 
Furthermore, if S4 7^ 0, DQ^d = 6) can be treated with a second BT iteration giving rise to logarithms. 



BaD, 



'o(d = 6) = J dS3V-'-'/\xuX2,xs) = J dS3[^ lny(xi,X2,X3) + l 



3 
1=0 



^ . 



(139) 



Note that each three-point function will have, after the second iteration, its own sub-leading BT coeffi- 
cient. As long as -B4 7^ 0, the case where some or all of the sub- leading B are zero is fully covered by the 
results of Section ^. Similarly to Eg. (|36[) , the case of G4 = det(ff) = reduces Dq to a combination of 
four Co functions. 



6.2 Form factors in the D-family 

Also form factors can be written in a compact way by using the functions 



2^ni,n2,n3= j dSg 4' 4' (^1 > ^2 , X3) , 

Ci(/jTl) = j dS2XiV-'~'l\[T+~'^). 



(140) 



Here we give one explicit example, 



^4 2^100 = \xiDQ{d = &)- Vwoid = 6) + i [(1 - Xi) Co(Ol) + ^ (X^ - X^+i) CiiiiTl] 

1=1 



where the d = 6 component gives 



(141) 



B4^Vioo{d = 6) = I dS3 (2x1 - -Xi) lny(xi,X2,X3) 



dS2 



(l-Xi) lnV{01) + Y,{Xi-X^+l)xl lnV{iT+i: 



i=l 



(142) 



Suppose now that i?4 w 0: an alternative derivation for Dq is needed. However, we first consider an 
extension of Eq.( |139| ). 
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6.3 Do at O (e) 

A three-point function with the insertion of a counterterm, which is needed in a two-loop calculation, 
is formally identical to a one-loop four-point function with one zero external momentum. Its expression 
at O (e) is needed and this follows from the BT algorithm quite nicely, 



We obtain 



IGBiD^^^ = 1- I dS3l3+ I dS2 



B, ^ 



B 



t=w,s,e,n 



3i 



i=w,s,e^n j=r,g,b 



B 



3i 



where the integrand is given by 

Jg = 2-h22 lny(xi,X2,X3) -h3 In^ X2, Xg), 



(143) 



dSiIuj, (144) 



(145) 



Ii = - Yl di\6lnVi{xuX2)+ln^Vi{xi,X2] 



i=w,s,e,n 



(146) 



/2I = 1 + 4 lnVi{xi,X2) + In^ Viixi,X2), luj = 2 Inl^j(xi) + In^ V^j{xi), (147) 
where we introduced the quadratics 

Vu, = V{xi,xi,X2), = y(xi,X2,x2), Ve = V {xi, X2, 0) , K = ^(1, 2;i , 2:2) , (148) 

= y,(xi,xi), V,g = Vi{xi,0), Vib = V,{l,x^), (149) 

and the auxiliary factors 

d^ = Xi-X2, ds = X2-X3, de=X3, d„ = l-Xi, (150) 

dir = Xii — Xi2, dig = Xi2 dih = 1 — Xii. (151) 
In the last equation Xi and Xij are the BT co-factors of the corresponding quadratics defined by 



Eqs.([l48D-([l49D 



6.4 Evaluation of Dq when ^ and Dq is regular 

This is the case when — > but the condition < X3 < X2 < Xi < 1 is not fulfilled, i.e. the point 
of coordinates Xi is outside the original simplex. Thus Dq is regular and the integration hyper-contour 
can be distorted to avoid the zeros of V without pinches. Two methods will be introduced to cover this 
case. 
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6.4.1 Method I 

We proceed analogously to the case of Cq. First of all we define Vb(xi, X2, X3) = V{xi,X2,X3) — B4 and 
write down the BT relations corresponding to -B4 7^ and i?4 = 0: 



1 + {X - X) X2, X3) = ^4 V-^-'/^{xi,X2, Xs), 



1 + 



1 



{x - X) 



Vr 



-l-e/2 







Xi,X2,X3) = 0. 



(152) 



2 + e 

Then we sum each side of the two equations and integrate by parts, obtaining: 

3 

Do = -l I dS3B^^V-^-'/\xi,X2,X3) I dS2Y.{X,-Xi+i)Bl^V-^-'/\fT+l) U, 

(153) 

where T^'^jsub = — ^o- The twofold integrals are nothing but differences of three-point functions 
that can be easily treated according to the results obtained for Cq in Section ^. For the first integral of 
Eq.( |153| ) we set e = and apply again the BT algorithms: 

l-^{x - X)da: lny(xi,X2,X3) = B4V''^{xi,X2,X3), 



Ba 



l + -ix-X)dx ^(xi,X2,X3) = 0, 

l-^{x - X)dx lnVo(2;i,X2,X3) = 0. 
By summing the first two equations and by subtracting the last one, we get: 



B^^V ^{xi,X2,X3) \snh= -^{X - X)d:^ 



lny(xi,a:2,a^3) 

-D4 



where we introduced the '++ distribution' 

'lny(2;i,X2,a;3)^ ^ In y(xi, X2, 3:3) - In V o(xi, X2, X3) - B4Vq'^{xi,X2,xs) 

++ 



R2 
^4 



R2 
^4 



Now another integration by parts is performed and we obtain: 



dS^B^^ V {xi,X2,X3) |sub= 2 



lny(a;i,3;2,a:3) 

-D4 



(154) 



(155) 



(156) 



(157) 



dS2Y,{X^-Xi+i) 



lnV{i i + 



i=Q 



-D4 



(158) 



In the twofold integral the '++' distribution contains terms such as V^~^(i i + which are three-point 
functions with e = and can be calculated in the usual way. In the first integral we have V(^^{xi,X2,X3) 
which can be reduced again to a sum of three-point functions: 



3 

/ dS3 Vq~\xi,X2,X3) = I dS2Y, i^i - ^i+l) V'(^'^l)- 



(159) 



This result is obtained by using the second relation of Eq.( [L54[ ) together with an integration by parts. 
An analogous procedure can be derived for tensor integrals. 
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6.4.2 Method II 



Also for the four-point function another method is available. First of all we perform a Taylor expansion 
around i?4 = in Eq.( |135D : 



dSsiX) 



-2-e/2 



Y,^n + l)V{n + 2){-B,T, 



n=0 



and then we use the BT relation 



1 + 



2n + e 



{x' H xy-'/^ = 0. 



(160) 



(161) 



In this way, an arbitrary coefficient in the Taylor expansion in i?4 can be reduced to a twofold integral 
to be treated according to a standard BT algorithm. There will be several new quadratic forms to be 
introduced, in two and one variables. We obtain the following form for the coefficients 



V{n)= f dS-i{X){x'Hx)-''-''^ = —^y XI: I dS2Q-''-''\Tii + l 

where X^ = Xi — Xj+i and the quadratic forms Q{Ti i + 1) are 

Q(rri) = Qo = g(l -Xi,xi- X2, X2 - Xs), 
Q{T r2) = Qi= Q{xi -Xi,xi- X2,X2 - X^), 
Q(r23) = Q2 = Q{xi - Xi,X2 - X2,X2 - Xs), 
Q(r34) = Q3 = Q{xi -Xi,X2- X2, -X3). 



(162) 



(163) 



In general we will write Qi = x^ Hi x + 2 Kj x + Li and introduce corresponding sub-leading BT factors 
Bi = Li — -f^j~^ Ki and co-factors Xi = —Kj with element Xij (j = 1, 2). A second iteration will 
bring into the result quadratic forms in one variable. 



Qi,o(a;) = Qi(l,x), Qi^i{x) = Qi{x,x), Qi,2ix) = Qiix,0), 



(164) 



which are written as Qij{x) = hij x'^ + 2kijX + kj and also have a sub-sub-leading BT factor h 
lij — kfj/hij and a co-factorXjj . By applying again the BT procedure we finally obtain for 2?(2) 



while for n ^ 2 

V{n) = 



V(2) 



1 + 



To JTo 



Po,ii(2)+ / dxVi,,j{2) 



InQij(l) 



InQy(O) - 2, Pi,y(2) = -\nQ,j{x) 



(3 - 2n) (n 



^ n ^ [^o.*(™)+ / dxVi,i{n)+ j dS2V2,i{n) 



1 



4(n-2) ^ 



i=0 

2 

E 

j=0 



hij 



-n+2( 



III J 



Ji] — 

(165) 
(166) 

(167) 
(168) 
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^^2(n-2)_^2n-5 



P2,,,(n) = ^Q-"+2(xi,X2). (169) 



Bi 



Now we can sum the series according to well-known results: 



X 



E 

n=l 

and obtain: 



where 



1 — X 



n=l 



^ (2n+ l)n 

n=l 



/" ln(l-xy2), (170) 

JO 



X, 



Do = V{2) +Y.^ + I dxDi^i + I dS2 D2,^ 



Do, 



2 x^- fi 



E 

j=0 



-'tj JO 



dy 



1 + 



(171) 



Q,,(l) + ^4y' - In Q^j(0)+-^4y' 



Qij(i) 



Q*j(0) 



I)2,*(n) 



L 

j=0 



1 



1 1 \ Q,Ax) + B^ ( 1 



Qij (2^) 



Qjj(x) +Biy^ 

Qij {x) 



2Bi 



Ba 



Qi{xi,X2) + B4 2 Qi(xi,X2) 



3 Qi(xi,X2) + -B4 ^ 3 ^^Qi{xi,X2) + B4,y'^ 



Qi{xi,X2) 



(172) 



In the twofold integral we still have a polynomial with a negative power; it is possible to increment the 
power of the latter by employing the BT relation: 



{Bi + S4) iQ^ + Bi)-^ = 1 - i (x - Xi) 5^ ln(Qi + B^) 



The usual integration by parts gives: 

1 



J dS2{Q^+Biy' 



2 {Bi + S4) 



2 

1 + 2 / dS2 HQi + B^)-Y, ^ [ dSi ln{Qij + B^) 
J j=0 J 



(173) 



• (174) 



In this way all coefficients in Eq.( |160D are explicitly computed. We can treat form factors with the same 
procedure; indeed the auxiliary functions 



3 

D{{n}) = / dS^\{x'I^V-^-^/\xi,X2,x^), 
•' 1=1 



(175) 



become linear combinations of 
P(n, 711,712,713) = 



dSo 



3 - 2 77 + Til + 712 + 713 
X {(Xi - X2r (^2 - (1 - ^1)"^+' Qo "(^1. ^2) 

+ (xi - xi)"i (xi - X2r {X2 - xsr {Xi - X2) Qr"(^i, ^2) 

+ {X2 - X2T' (X2 - X^r^ (Xi - Xi)"i {X2 - X3) Q2 "(^1, 3^2) 
- (Xi -Xi)"i (X2 -^2)"^ (-X3)"^+1Q3"(X1,X2)}, 

which are form factors of the C-family. 



(176) 
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6.4.3 Sub-leading behavior 

The above derivation assumes that the sub-leading BT factors are not zero. If this is the case we start 
from the fact that Dq is a sum of terms of the form 

°° -I- 1 r 
^0 = - E (- ^4)" / dS2 Q-"-^ (177) 

n=0 

with Qi constant. When ~ we cannot apply the BT algorithm and it is necessary to adopt an 
alternative strategy. If the generalized three-point function associated to the vanishing coefficient Bi is 
regular, we can perform another Taylor expansion exactly as it was done for Cq. Otherwise we rewrite 
the sum as follows: 

E ^ Q--' {-B.T = {B.)-'" dz (Q, + z')-\ (178) 

n=0 •'^ 

To the price of introducing an extra integration we are now in the situation where the series is replaced 
by a single term and the sub-leading BT factors for the xi — X2 integration are replaced by Bi + and 
the singular behavior is left in the last z integration. Note that, if -B4 < we simply have to replace 
everywhere B^ and Qi with — i?4, —Qi. Actually the best way of dealing with the singular behavior is to 
change variable and obtain 



j dS2{X) dz{x'HiX + B3, + B4Z^)-^, 



(179) 



and perform a Mellin-Barnes splitting (accompanied by a sector decomposition) of Hi x and B^i+B^ z^ . 
The leading behavior comes from the single pole at s = 1 and it is, therefore, proportional to 



(S4i?3.)-'/'arctan(-^) , (180) 
where B ^ B — iS. Using the identity of Eg. ( [1 781 ) in Eq. (|177] ) we obtain 

D'o = -q^ C dz [ dS2{Qi + B4Z^)-\ (181) 



JO J 

The first step toward the evaluation of the integral consists in splitting the integral in the sum of two 
parts: 

Di = l{-li li;2 = -Qi f dz f dxi [ dx2iQi + B4zY^. (182) 

Jo Jo J0;xi 

Now the first integral is singular when B^i while the second is not. Next, we consider the calculation 
of /|. Introducing the new variables of integration x'^ = Xi — Xi, we have that 



If = -Qi J dC2{X) dz{x'Hx + B3, + B4 z2)-2 



(183) 



then we apply the Mellin-Barnes transform and obtain 

»1 fioo 



^l = -Y-i I ^^^W / dz I dsB{s,2-s){x'Hx)-'pl^^, (184) 
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where l/yCmix = B^z'^ + B^i. It is possible to perform the integration over the variables Xj and the result 
can be cast into the following, general, form: 

f{s) 



J dC2{X){x*Hx) 



s-l 



As a consequence we obtain 



2m 



1 ^ioo 

dz 



dsB{s,2-s)p^^^M. 



(185) 



(186) 



We can now perform the integration over s and the leading divergent behavior is given by the pole at 
s = l: 



Ij'-°g = g, /(l)arctan 



-1/2 /S4^ 



B 



3i 



(187) 



where the determination of the arctan function is fixed by the —i5 prescription. The regular part of the 
integral /{ is 



00 „i 

/j'-g = -q, V (-1)" / dz f{n + 2) {B,z^ + i?3i)" , 

n=0 

the integration over z can then be performed: 

00 ^ 



n=0 



2 ' ' 2 ' 5 s 



3 5^ 
hi'' 

The integral /2 can be calculated by performing a Taylor expansion around i?42;^ + -634 = 0: 

00 1 1 „i 

4 = -ftE / dzi-B^z^ -Bsi)'' dxi / dx2(fe + l)[(x-X)*i//- 
i._n -^0 •''0 Jxi 



(188) 



(189) 



k=0 

00 



1 3 R 



fc=0 



where we introduced the integral T2{n) defined by 

l2{n) = [ dxi [ dx2 [{x - XY Hi (x - X)]-" = / dS2{X) 

Jo Jxi 



X* Hj X 



(190) 



(191) 



These integrals can be treated in analogy to what was done for C{n) (Eq.(|54)) and D(n) (Eq.( p.62]) ); using 
the relation 

(x*/?x)-" = 0, (192) 



1 + 



2n 



we obtain 



2{n 



^ 1^' dx [(1 - X2) + {X2 - X,) r^^ + X, r^^] , 



(193) 



where P(xi, X2) = x^ H x and where 

Vi=Vil-X2,x-Xi), V2 = V{x-X2,x-Xi), V3 = Vix-X2,Xi). (194) 

The remaining one-dimensional integrals are now generalized Bq functions which can be calculated em- 
ploying the methods introduced in Section |^. 
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6.5 Evaluation of Dq when B4 w and Dq is singular 

This is the case when ^4 — > and the condition < < X2 < Xi < 1 is fulfilled, i.e. the point 
of coordinates Xi is inside the original simplex. Thus Dq is singular and the integration hyper-contour 
cannot be distorted due to the occurrence of a pinch. In order to minimize the number of terms in the 
solution we write 

1 3 »i 

Yl dxi = 2_\ I / / dxk- (195) 

1 Jo Jo Jo 

1=1 cycl.pcrm. 

Furthermore, for each term we may use the following identity 

XI pxi pxi pX2 rxi px^ 

dx2 / dxs = / dx2 / dxs + / dx^ / dx2- (196) 
'0 Jo Jo Jo Jo Jo 

Therefore the integral over [0, 1]^ can be written as 

r H dx,V-^~'/\xi,X2,xs) = [ dS^Y. y~^~"'^{xi,X2,x^), (197) 

•''^ j,=l perm 

where the sum is over all permutations of {xi,X2,X3}. The first term in the sum is exactly our original 
Dq function while the rest gives the five complementary functions which, by construction, are regular 
and can be treated according to the strategy of Section 6.4. As a consequence, we now have to evaluate 
D'^^'^ when i34 ~ and the point of coordinates Xi is inside the unit cube. Using standard techniques, 
already applied to the Co-function, we obtain 



i=l J~^i 

dxi I dx2 I dxs {Qi + B^y^-'/^ (198) 
.■1 Jo Jo Jo 



i=l 



If (3(xi,X2,X3) = X H X the quadrics Qi are defined as 

Ql = (5(xi,X2,X3), Q2= Q{-Xl,X2,X^), Qs = Q{-Xl, -X2, X3) , 

Q4 = Q{-Xi,-X2,-X3), Q5 = Q{-Xl,X2,-X3), Qg = Q{xi, -X2, X3) , 

= Q(xi, -X2, -X3), Qs = Q{X1,X2,-X3), (199) 

and moreover the coefficients , /5j and ji are defined by 

Qi = ag = 07 = Qg = 1 — Xi, a2 = 03 = 04 = 0:5 = Xi, 

/3l = /?2 = /?5 = /?8 = 1 - X2, P3 = P^ = Pq = I3j = X2, 

71 = 72 = 73 = 76 = 1 - -^3, 74 = 75 = 77 = 78 = -^3- (200) 
In this way the evaluation of 1)™^° is reduced to the one of 

r^x rPi pi 

Di= / dxi / dx2 / dx3{Qi + Bi)-'^-'/'^, (201) 
Jo Jo Jo 
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for which we introduce a MeUin-Barnes transform 

Di = 7r- dspl-'Viis), (202) 

where /34 = l/i?4. It follows that 

T>i{s) = aiPijiB{s , 2- s)di{s), di{s) = J dC3Q~^{aiXi,/3iX2,jiX3). (203) 
After a sector decomposition in [0, 1]^ we obtain 



dC. 



Qj_''{ai,PiX2,-yiX3) + X2 Qi''{aiXi,Pi,-jiX3) + x^ Qi''{aiXi,/3iX2,'ji) 



(204) 



Collecting all pieces we obtain 

8 3 

i=l j=l 



27ri 



+ j oo 



J dC2Q-^(,xi,X2). (205) 



The evaluation of D^^^^ has been therefore reduced to the one of an inverse Mellin-Barnes transform: 



D 



3-2s 



(206) 



For large values of \p\ we close the s-contour over the right-hand complex half-plane at infinity with a 
single pole at s = 3/2 and also ats = 2-|-A;,A;>0. The result is straightforward 



4 ■' "-2' 



k=0 



2k + l 



i-p)- 
k\ 



Since / stands generically for 



/(S)~ J dC2Q-j'{x,,X2), 



(207) 



(208) 



we have to evaluate this integral for s = 3/2 as well as for s equal to an arbitrary integer > 2. We will 
use the following two recursion relations: 

J dSiQ-\xi) = j^-^j^[{2s-3) J d5iQ^-^(xi)+XiQi-^(0) + (l-Xi)Qi-^(l) 



J dC2Q-'{xi,X2) = ^j^^jy^{2{s-2) J dC2Q'-'{xi,X2) + J dSi [x2Ql-%xi) 
+ Xi Ql-'ixi) + (1 - X2) Ql-'ixi) + (1 - Xi) Ql'^xi)]], 



(209) 



where Q is an arbitrary quadratic form in one (two) variables, B is the associated BT factor, X the 
associated BT co-factor and Qi are secondary quadratic forms defined by 



Qi{x) = Q{x,0), Q2{x) = Q{0,x), Qsix) = Q{x,l), Qi{x) = Q{l,x). 



(210) 
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We give in the following the explicit solution for s = 3/2. 



/ 



dC2Q 



-3/2, 



I J dC2Q'/\xuX2) + J d5i/i(^)-/o(^)], (211) 



where we have 



/i(^) = ^2 ( + ^) (xi) + ^1 (;| + ^) (^i) 



+ (1 - X2) (J- + ^) Qf (xi) + (1 - Xi) (J- + ^ 



X 



/o(^) 



^iiQ}^'(0) + (l-Xn)Qi'"(l 
1 -X2 



+ 
+ 



+ 



Bo 



(xi), 

X12 Ql^\0) + {l-Xi2)Ql^\l] 



B3 
1-Xi 



^13Q3^'(0) + (l-Xi3)Q3/'(l) 

gy'(o) + (i-Xi4)gy'(i) 



(212) 



where B and Xi , X2 are the leading BT factor and co- factors of the original quadric while Bi and Xu are 
the sub-leading ones for the secondary quadratic forms. Note that in this case the BT procedure allows 
us to derive a result in terms of roots of quadratic forms. For s = 2 we derive instead 



/ 



dC2Q '^{xi,X2) 



1 

Jb 



d5i/i(2)-/o(2) 



(213) 



where we have 



/i(2) = ^ InQi(xi) + ^ lnQ2(xi) + InQ^ixi) + lnQ4(xi), 

±)1 ±)2 ^3 -D4 



/o(2) 



+ 
+ 



X2 
Bi L 
I-X2 r 

B3 
l-Xi r 



+ 



Ba 



Xn InQi(O) + (1 - Xn) InQi(l) - 2 
Xi3 lnQ3(0) + (1 - X13) lnQ3(l) - 2 
Xi4 lnQ4(0) + (1 - X14) lnQ4(l) - 2 



Xi 



B' 



2 L 



X12 lnQ2(0) + (1 - X12) lnQ2(l) - 2 



(214) 



For s > 3 we derive similar expressions, not to be reported here. Dq'^^'^ is, therefore, explicitly known. 

The procedure cannot work if one of the sub-leading B is approaching zero. Again, the B are sub- 
leading of second type as explained in the discussion before Eq.(85). In the case i? ~ the solution consists 
in a double Mellin-Barnes splitting. Consider Eqs.(206)-(208), the evaluation of D™^'^ is reduced to the 
one of 



D 



1 



2iTi 



+i 00 



B{s,2-s) 
3-2s 



dC2Q '{XI,X2), 



(215) 



where Q denotes one of the quadratic forms of indices ij, which we write as Q{xi , X2) = x^ H x+2 x+L. 
Let BsL = L — K, psL = ^/Bsl {Pl = Pi) and X = —H~^ K. After transforming xi,X2 we obtain 



D 



27ri 



+i 00 



,_, B{s^2-s) 
'^"^ 3 -2s 



dC2{X) {x^Hx + Bs 



(216) 
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With < Re t < Re s < 2 we write 



D 



+ 1 oo 



, , 9 „ „ / B{s, 2 — s) 

ds dt «r — - 



B{t,s-t) J dC2{X) {x^H: 



(217) 



The xi,X2 integral has been already considered, starting from Eq.(^), and we know that it develops a 
pole at t = 1 times the integral of a quadratic form in one variable. If, for the latter, the BT factor is 
not around zero we can safely proceed. If we write 



dC2{X) {x^ Hx) 



-t 



fit) 



then we get 



D 



27ri 



+i oo 



+ pI-' 



L psL 3 -2s \ Pl 



k=0 



3-2s 



Next we use k = Psl/Pl and 



27ri 



+j oo 



ds 



r(2-s) r(s- 1) 

3 - 2s 



E 

k=0 



oo 

K \ - 

9 2^ 



-1)'' T{k + l) 



kl 



kl 2k + l 



k=0 



r (fc + 1) r (fc + 1/2) 

r(A: + 3/2) 



kl 



K^/^ arctan(K" 



-1/2N 



Therefore we get 



(218) 



(219) 



(220) 



D 



(p, ps.f" /(I) arctan [P^]"" + \ p^fl' E ^ ( ^ + ^ ) + ^ 



PSL 



k=0 



-PSL 



A:! 



E 

k,l=0 



r{k + i + i) f{k + i + 2) (-psl)-'' (-pl)-' 



21 + 1 



A;! 



II 



(2211 



In the above results it is always understood that Pl ^ Pl — and that pgi^ — > psL — i6. 

As illustrated in the previous section the evaluation of f{k) is reduced to computing a one dimensional 
integral of a quadratic form in one variable (with a small negative imaginary part) which can be done by 
BT iteration, unless the corresponding BT factor is around zero. If this is the case, we start wondering 
why the world is against us or return to the starting expression and use one more Mellin-Barnes splitting. 



6.6 A new integral representation for Do 

Also for Dq we are able to present a new integral representation which is suitable for numerical 
treatment, especially when = 0. Consider again Eq.(201): after a sector decomposition in [0, 1]^ we 
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obtain 

Di = ai Pi 7i 



xl Qi{ai, Pi X2, 7i X3) + B4-i6 



-2-e/2 



+xl 



xl Qi{aiXi,Pi,jiX3) +B4-iS 
3 



-2-e/2 



+ xi 



dCs X3 Qij {xi,X2) + B4^ — id 



xl Qi{ai xi,Pi X2, 7i) +B4-i6 

2-e/2 



-2-e/2 



(222) 



where we introduced the quadratic form 



They are expressible as 



Qii{xi,X2) = Qi{aiXi,l3iX2,^i), 

Qi2{xi,X2) = Qi{aiXi,l3i,JiX2), 
Qi3{xi,X2) = Qi{ai, Pi Xi,%X2). 



Qij {xi , X2) = X* Hij X + 2 Kfj X + Lij = Mij x, 



where = {xi,X2), = {xi,X2, 1) and the matrix M is given by 



Li, 



Now we apply the results for the master integral of Appendix A with coefficients 

a = Qij{xi,X2), h = Bi, a = 2 + e/2, /3 = 2. 
With the usual replacement Qij, B4 — > Qij, B4 — i 5 we obtain 



(223) 

(224) 

(225) 
(226) 



Di = a, 



P^1^Y.J 1 1 B-'/'-'/' i? 0, i + £ j Q7f/\x,,X2) - dx^ X^ Q-;^-''\x^, X2, X3) 



(227) 



In the last equation polynomials Qij{xi,X2) and (5ij(xi, X2, X3) = Qij{xi,X2) + -B4 x\ make their appear- 
ance; their powers can be easily incremented with the BT method, as long as the corresponding factors 
B^ij are different from zero. The BT relations to be used for this purpose are as follows: 



Qij ^{xi,X2) 



B: 



I 2 

1 - 2^ ^{Xl-Xiij)d:ci 



1=1 



Qij{xi,X2), 



Qfj ^{xi,X2,X3) 



1 



B 



3ij 



1 / 



Qij{xi,X2,X3), 



where the sub-leading BT-factors are defined by 

det(M,,) afpfjf 



B3 



G: 



3ij 



G: 



3ij 



det{H), Guj = AssiMij), 



(228) 
(229) 

(230) 
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X 



lij 



Ai3(Mij) 



A23(Mij) 



Mij - — , [2.61) 

where Aj^.(Mjj) is the co-determinant of the element Ik of the matrix Mij. The BT method cannot be 
used when also some of the B^ij are near zero and, in this case, we have to proceed along the same lines 
used for the three-point function. We rewrite the polynomials Qijixi, X2) in the following way: 



Qij {xi,X2) = (x - XijY Hij (x - Xij) + Bsij. 



(232) 



Next we transform variables according to xi = x\ + Xnj , x^ = x'^ + and split the x\ — X2 integration 
interval according to 

dxi=(- + dxi. (233) 

-XUA ^ JO ' 



Each of the integrals is mapped into the unit square and we derive 

3 4 



3 1 e 
2'2 2 



Qijl{xi,X2) + B^ij 



-3/2 



Jo 



dX3 x| 



Qijl{xi,X2) + Bsij + B4 x| 



-2-e/2 



(234) 



where the replacement B^ij B^ij —i6 is understood and where we introduced new auxiliary quantities: 

llijl = Vij2 = 1 — Xiij, 1]ij3 = r]ij4 = Xiij, 

Cijl = Vij3 = 1 — X2ij, iij2 = iiji = ^2ij- (235) 

The new quadrics are defined by 

Q{ij)l = Q'ijii^ - ^lij) Xl, (1- X2ij) X2), Q(ij)2 = Q'ij((l - ^lij) Xl , -X2ij X2), 

Q{ij)3 = Q'iji-^Uj Xl, {I- X2ij) X2), Q(ij)4 = Qiji-^lij Xl , -X2ij X2), (236) 

with Q'-j = X* Hij X. At this point, after performing the usual sector decomposition of the square [0, 1]^, 
we apply again the master formula of Appendix A and get: 

3 4 2 „i 

(237) 



d 1 1 
A = ttj A7i X] X] ^ rjiji^iji / dxiDijijn, 

j=l 1=1 m=l 



Dijim "2-^4 ^31^^ ( 2' 2 2 ) ^v^i^^^^^ 



v2' 2 2 
-^i?(l,l + e/2)Q-j„(xi)^'dx3X 



B4 x| + Bsij 



-l-e/2 



1 1/2 e/2 ^ Q 1 ^ ^ QJ^{xi,X2) + dX2 ^ dxg x| Q J^'/^(xi, X2, X3), 

(238) 



where we introduced Qijiiixi) = Qiji{xi, 1) and Qiji2{xi) = Qiji{l,xi). Likewise we have 

Qijlm{xi, X2) = Qijlm{xi) + Bsij xl, 
Qijlmixi,X2,Xs) = Qijlm{xi) + B^ij x| + B4 x|. 



(239) 
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The divergent behavior of the first and of the third term is evident and the second integral depends 
essentially on the ratio 

- ^'-'^ (240) 



Therefore, by applying the master formula of Appendix A, we obtain 



-B4 X2 + B-sij 



-l-e/2 1 



1/2 „-l/2-e/2 



BZ7 B 



B['-+'-,'- 
'2 2' 2 







/ dx3 




Jo 





2 3ij 4 

I i?3,]/' Bl'" - B,' arctan(K,r.V2) + O (e) 



-l-e/2 



(241) 



Note that the first term in Eq. ( |24lD cancels the first term of Eq. ( P38D giving our integral representation 
for A: 

3 4 2 1 j-^ 

A = A 7i X] ^^i' ^^i' / 2 ^^^^^ ^ arctan(K,r>/2) QT.]^{xi) 

j=l 1=1 m=l I 

-jB'^^'^J^ dx2Q~/il^{xi,X2) + / c^a:3Q^:;^|^(xi,X2,X3)|. 

(242) 

As a final step, the power of the polynomials appearing in Eq.( ^42D can be incremented by using the 
following set of BT relations: 



1 



^ijlm 



(243) 



Q^jlL{xi,X2) 



1 f 1 

~ 271 



-'ijlm 



xi + ) + X2 



ijlm 



Qiilmixi,X2) 



2/i 



Hjlm 



fc=2 



where we defined 



Qijlmixi^ hijijYi 2^1 ~l~ 2 k^jifYi Xi -\- lijlm ^ijlm ^ 



Qf,7m(^l'^2,X3), 



A-2 



ijlm 



i]lm 



h 



ijlm 



If also bijim tends to 0, we must repeat the procedure for the quadratic form Qijim{xi) 



(244) 
(245) 



(246) 



7 Infrared divergent Dq 

At the one-loop level the evaluation of infrared divergent A^-point functions with A^ > 4 is greatly 
simplified insofar the divergent parts can always be reduced to three-point functions. This decomposition 
is even more relevant in our approach where the evaluation of the infrared finite remainder does not 
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require any additional work. Let us give an example for Dq [^], where we consider the configuration 
with 



2 2 
Pi =P2 



2 2 2 

-"^2, Ps = Pa 



-771,4, '^1 ~ "^3 ~ 0- 



If we introduce the following notation: 

[l]=q\ [2] = {q+pif + ml [3] = {q + Pf, [i] = - p^f + rnj, 
where P = pi + p2, then the following decomposition holds 



7d ^ ' il + P) ^ // , n2 2 2 n ^ 

da - Co{{p2+P3) ,P2,P3; m2, 0,1714) 



[1] [2] [3] [4] 



- Co{pI,pI, {P2+P3f ; 0,m2,m4), 
showing two IR divergent Co functions and an IR finite integral. For the latter we obtain 

1 



T-jrem 

^0 ~ ~ 9 



j dS^V-^-'l^ {k^ - k-P + 2V), 



where 



and the quadric 



[1] [2] [3] [4] 

3 

P = Pl+P2, Q=P2+P3, k = ^XiPi 



1=1 



V = ml {xj + X2) +m\xl- (P^ + 2 mf) xi X2 + (Q^ + P'^ +ml + m\) xi X3 
- {Q^ + ml + ml) X2 X3 + {x2 - X3), 

shows a zero for xi = X2 = x^ = 0. The numerator in Eq.( p50D becomes 

1 



k"^ - k ■ P + 2V = V - - P^ {xi- X2 + X3), 



so that we obtain 



7-) rem 
^0 



dS3 



V-'-^'^-\p^ {X,-X2+X3)V-^-^/^ 



(247) 
(248) 



(249) 



(250) 



(251) 



(252) 



(253) 



(254) 



By infrared power counting the expression for Dq*^™ is manifestly finite and can be computed according 
to standard BT techniques when 7^ 0. The explicit expression of the coefficient -B4 is 



1 P^ 



dei{H) 16 



4 m\ m\ — (Q^ + m2 + 771,4) 



2\2 



1 P 



\- A(-(3^, 77x2,7714), 



det(i?) 16 



where det(i?) is the determinant of the matrix H in the BT relation: 



det(iJ) = — [A(-g2, ml,m\) + 



(255) 



(256) 



As a consequence P4 = only for < 0, i.e. for normal(pseudo)-threshold = —(772,2 i 7714)^. 
Around the thresholds other techniques are used since the BT procedure overestimates the singularity, 
see Section 6.4 for the pseudo-threshold (where Z)q^™ is regular) and Section 6.5 for normal threshold 
(where Dq^'^ is indeed singular). 
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7.1 Classification of IR divergent cases for the D-family 

A short classification of IR divergent four-point functions can be based on IR power counting which, 
for one-loop diagrams, requires f^d'^q/q'^. The configurations of interest are those where some of the 
internal masses vanish and the Mandelstam invariants of the process are not constrained. 

Without loss of generality we may assume that mi = 0. For q = A and A vanishing the requested 
behavior of O (A^) is obtained in three different cases Q 

1. pI = -ml pI = -ml; 

2. pI = -ml = — m|; 

3. pI = -ml p2 = -ml- 

Note that only the first one is of interest, so the IR configuration is 

mi = 0, pI = -ml, pI = -ml- (257) 

If also ms = then there is an additional pole for q = —P and we introduce q = —P + A. The four 
denominators in Dq behave as 

[1]-P2^ [2] ~pi + m2-2p2-A, [3]-A2, [4] ~ + ml + 2p3 • A, (258) 

and the requested behavior of O (A^) is obtained ii pi = —ml and p| = —ml, as expected. Therefore 
another IR configuration is 

mi = ma = 0, {p\ = —ml pi = —ml} and/or {pi = —ml p\ = —ml}. (259) 

If instead we have mi = and m2 = there will be another pole a,t q = —pi whose IR nature depends 
on p2 = Q and pi = —fn\- The new IR configuration is 

mi = m2 = 0, Pi = 0, {pI = —ml} and/or {pi = -m^}. (260) 

Furthermore, if mi = m2 = ms = there are the following IR configurations: 

mi = m2 = m3 = 0, {pj = Opl = -ml} or {pj = pi = 0} or {pi = pj = -ml}. (261) 

In all cases the IR box can be reduced to a combination of (possibly) IR Co functions and a IR finite 
reminder. 

To further clarify the connection between IR configurations and Landau singularities we consider the 
Dq specified by Eq.( |247D . The corresponding Landau equations are 

aiq'^ = 0, a2[{q+Pi)^ + ml]=0, 03(^ + ^)^ = 0, [{q - p^)'^ + ml] = 

ai q^l + 02 (g + pi)^ + 03 (g' + P)^ + 04 (g - p^)^ = . (262) 

The equations in the first row of Eq. ( |262| ) are satisfied, for Oj ^ 0, if 

^2 = 0, g-pi=0, q.p=-]^P\ q-p^ = 0. (263) 
^This is the well-known eikonal approximation [E2|, whose validity we assume without further justifications [E3|. 
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By contracting the free Lorentz index in the equation in the second row of Eg . (12621 ) by g^, (pi)^, P^, and 
P4^, and by using Eg. ( p63D we obtain the system of eguations 

as = , —2 02 + as — 2 04 • P4 = , 
-ai p2 + as = , 2 a2 pi ■ P4 - as P^ + 2 a4 ml = . (264) 

If the system of Eg .( 12641) admits a proper solution (i.e. a solution in which not all the aj vanish), the 
determinant of the system must be egual to zero. One can readily verify that the determinant is 



A 



16 



\{-Q^,ml,ml). 



(265) 



Therefore the leading Landau singularity reguires < and corresponds to the usual normal threshold 
= — (m2 + mif' on top of the infrared pole. 



8 Five-point functions (E-family) 

Several results are known in the literature about the reduction of dimensionally regulated one-loop 
A^-point integrals with > 5 ([^, [24| and [^]). The generic, scalar, pentagon is given by 



where H, 



Eo = j dS4V-^"'/^{xi,X2,X3,X4), V{x) = x^Hx + 2K^x + L, 
- Pi ■ pj with z, J = 1, . . . , 4, L = ml, and where 



Ki = 


1 

~2 


{ml 


- ml 


-Pi -Pi), 






K2 = 


1 

~2 


{ml 


-ml 


- 2pi ■P2 


-P2 ■P2), 




K3 = 


1 

~2 


{ml 


- ml 


- 2pi -ps 


- 2p2 ■P3 


-P3-P3), 


K4 = 


1 

~2 


{ml 


-ml 


- 2pi -^4 


- 2p2 -PA 


- 2ps -pA-PA - Pa) 



(266) 



(267) 



mi 




Figure 4: The one-loop, five-point Green function of Eq.(266). Propagators are q -\-m\ ■ ■ ■ {q+pi 



\-pi)'^ + ml. 
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8.1 Evaluation of Eq when 7^ 

It is a virtue of the BT algorithm that we can show the decomposition of £'0 in five boxes (in d = 4) 
with just one iteration, as depicted in Fig. ^. We obtain 

^o = rB- E^^^o^' (268) 

where the weights are 

Wi = X,-Xi+i, Xo = l, X5 = 0, (269) 

and where = L — K and X = —K^ H^^. The further advantage in this derivation is that the 

nature of the weights is transparent since i^s = corresponds to a Landau singularity of the pentagon. 
Furthermore the boxes are specified by 

= j dS:,V'^~''^{ii + l), (270) 

where the contractions are 

1 = (l,xi,X2,X3), 4 5 = (xi, j;2,a;3,0), (271) 

etc. As long as is not around zero the derivation for the pentagon is completed since we know how to 
deal with boxes even if their (sub-leading) BT factors are around zero. 



8.2 Form factors in the E-family 

The nice feature of the scalar pentagon is that fourfold integrals disappear in the final answer, Eq.( |268| ). 
The same is not immediately true for form factors. En etc, which contain powers of the Feynman 
parameters in the numerator. We can use a new identity, 

{xi - Xi) V>'{xi,X2, X3, X4) = J^l^^ dj V>'+\xi,X2,X3,X4). (272) 

With the help of the identity and after integration-by-parts we are able to remove again all fourfold 
integrals in the form factors of the E'lj-series. To give an example we consider 

£^11 = J dSi XI V~^~'/\xi,X2,X3,Xi) = 4^ / dSs {E^^ - ^ Ef/). (273) 

Secondary quadrics are defined as: 

Va{xi,X2,X3) = y(xi,Xi,X2,X3), Vb(xi,X2,X3) = ^ (xi , X2 , ^2, X3) , 
Vc(xi,X2,a;3) = y(xi,X2,X3,X3), Vd{xi,X2,X'i) = y(xi,X2,X3,0), 

V;(xi,X2,X3) = y(l,Xi,X2,X3). (274) 

The Ell form factor is fully specified by 

E^^ = {Xi - X2) XI V~^-''^ + {X2 - Xs) XI V-^-'^^ + (X3 - X4) XI 

+ X4 XI V.^-'l^ + (1 - Xi) V-^-^'\ (275) 
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+ H^l V-'-^'^ - i?r/ V-'-^'\ (276) 
Similarly we obtain 

+ iiTl - V-'-^/\ (277) 

for i = 1, . . . , 4 and 

E^^^ = (Xi - X2) XI V-^-''^ + {X2 - X3) X2 y-'-'/' + (X3 - X4) X2 

+ X2 y-^-''^ + (1 - Xi) XI 

^^gl^ = (Xi - X2) X2 + (^2 - ^3) ^2 V-''-'''' + (X3 - X4) X3 

+ X4 ^3 F^'-'/' + (1 - Xi) ar2 

El^ = (Xi - X2) X3 y-'-^/' + (X2 - X3) X3 + (X3 - X4) X3 V-^-'/^ 

+ {l-Xi)x3V-^-'/^. (278) 

As it has been done before and for the reader's convenience we list all the form factors in the ii^-family. 
Note that there are also delicate points connected with the decomposition into objects belonging to the 
family since the decomposition itself is strictly defined in 4 dimensions. Furthermore, starting from six 
powers of momenta in the numerator, we will encounter UV divergent terms so that some care is needed. 
Clearly we can have such a term while working in a general ii^-gauge. 
With two momenta in the numerator we define form factors as 



E2;i 


= E2{ii), 




E2;5 


= E2{12), 


E2:6 


E2;9 


= ^2(24), 


E2;10 



^2(34), (279) 
where the auxiliary function E2{ij) is 

E2{ij) = J dS^XiXjV-^-'/\xi,X2,X3,Xi). (280) 

However, a new form factor arises, the one proportional to S/^^,, 

E2;ii = \ j dSiV-^-'l^{xi,X2,X3,x^). (281) 

£?2;ii can still be reduced to form factors of the £'-family by writing 

4 4 

= 4 [ Hij E2 {ij) + 2j2KiEii + L Eo\ . (282) 

i,j=l i=l 

Let us give the complete expressions for the form factors of the £?2-series. 

4i?5 E2{ij) = - j dS^ H;.' y-i-^/2(xi, X2, X3, X4) + j dSs {E^;^j - E^^.), (283) 
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where the various coefficients are: 



^2;ll 




-X2) 




'" + {X2-X,)xlV-'-''' 




+ (^3 




xfV-^- 


+ X4 X? 1^'-^/' + (1 - Xl) V-'-^l\ 


-^2;22 


= (Xi 


-X2) 




^i' + {X2-x,)xlv;'-''' 




+ (^3 






+ X, xl V,'-'l' + (1 - Xl) xl 


p-2 
-^2;33 


= (^1 


-X2) 




^/^ + (X2-X3)xiF-^-^/^ 




+ (^3 


-X4) 




-/^ + X, xl V-^-^l^ + (1 - Xl) xl V-^-^l\ 


-^2;44 


= {X, 


-X2) 




-^l^ + {X2-X,)xlV-^-^l^ 




+ {X3 


-X4) 


4v-'- 


-/2 + (l-Xi)ariF-2-e/2^ 


p-2 
-^2;12 


= ix. 


-^2) 


^1 


-/2 + (X2-X3)a;iar2V'-'/' 




+ {X3 


-X4) 


Xl X2 


-2-6/2 ^ js^^ V^^"'^^ + (1 - Xl) Xl 



(284) 

I XI y-^-/^ + H-^l XI v-'-^'' - vr'-^i\ 
I ^1 K-'-^/' + (^2-2' - H2i) X2 v;'-''^ 

I X2 V-^-^'^ + H^l X2 V~'-^l' - H-^l XI V-'-^l\ 

^ X2 V-'-^l' + (^32^ - H^i) X2 V-^-'l^ 

I X3 V-'-^l' + H^l xs V-'-^/' - H,,' X2 V-'-l\ 

' ^3 Va^-'l'' + (^4-2' - ^4-3') ^3 V'"^' 

^x,V-^-^l^-Hllx3V-^-^l^ 

I Xl + [(ij-1 - i/-/) X2 + {H22 - H2i) Xl] v;'-''^ 

{Hii - H^l) X2 + {H^i - H:^1)xi] 
+ [H^l X2 + Xl) y-^-'l'' - (H^,' Xl + H^,') V-'-'/\ (285) 
All functions Eij contain a term which is proportional to EQ{d = 8). Consider again E2-ii which becomes 

4^5 £;2;ii = -Eo{d = ^) + \ j [{Xl - X2) v-^-'/^ + {X2 - Xs) y-^-'i'' 
+ (X3 - X4) ^ v7'-'/' + (1 - Xl) v-^-'i^ 

When we consider 



-^2^11 = (-f^n^ - 

+ {H^i - Hil 

-^2;22 — {^21 ~ ^22 

+ (-f^23^ - -f^24^ 
-^2;33 ~ (-^3/ " ^32 

+ (-^33^ ~ -^34^ 
-^2;44 = (-^^41^ " ^42 

+ (^43^ - ^44^ 
2-^2;12 — i-^ll ~ ^22 

+ 



E, 



1 



iTT^ r (3) 



{q^ + ml) ■■■ {{q+pi+ ■■■ + Pa)"^ + 



2V 



(286) 
(287) 



the contribution proportional to £^0 {d = 8) is 

4 



Epjd = 8) 

4^5 



V+ 1^ Hi/Pii^Pji^ 



+ form factors N < 5. 



(288) 
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Since the four- vectors span d = A space-time the term proportional to Eo{d = 8) disappears and, 
therefore, there is no need to compute it. Indeed let us contract with pi, I = 1, ... ,4 to obtain 

4 4 

Piu {Sf,. + Yl Hr.^Pii^Pj-) = Pii^ - Yl ^ij" Pii^ = 0- (289) 

i,j=l ij=l 

When we go to three powers of q in the numerator the following result is valid: 
E^ua = - J dSi xi {5pi}f,^,a + Y {PiP3Vl}i^^oc + form factors N <b. (290) 

where the coefficients are 

elji = Hr.^xi + cyclic, (291) 

and where we introduced 

{Spj^iua = Pa + Sfj,a Pv + Pfi, {pqk}fj,ua = Pfj, Qu K + Cyclic, (292) 

showing that EQ{d = 8) disappears again. It is straightforward to extend the demonstration up to five 
powers of loop momentum in the numerator. This is the maximum of non-contracted powers that we can 
have in a renormalizable theory and starting from six powers we will have numerators of the following 
structure, 

•••'?A'5> (293) 

where scalar products can be simplified according to = [1] — m\ etc. There is only one case where 
the argument fails: suppose that one external line of momentum pi splits into two lines of momenta 
Pij,j = 1, 2. Then scalar products q ■ pij do not occur in any of the propagators and in the final answer 
we would end up with evanescent operators like 



i<5p^j. (294) 



We now describe the solution around = 0. 



8.3 Evaluation of Eq when f» and Eq is regular 

If S5 and the condition < < Xj+i <l,i = l,--,4is not satisfied Eq is regular but the 
decomposition into a sum of boxes fails. We write 

Eo = J dS4iX) {x^Hx + 55)-3-^/2. (295) 

Since Eq is regular at -B5 = we perform a Taylor expansion in with coefficients £{n), 

^ 00 

£^0 = 2 E + 1) + 2) ^(^ + 3) (- ^5)", (296) 

n=0 



with n > 3. Using 



£{n) = J dS4{X){x^Hx)-''-^/^, (297) 
J dSi{X) [1 + 2^^] {x' H xy-'/^ = 0, (298) 
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we easily obtain 



1 ^ r ^ 

f(n) = j32-E / dSs{Xi-Xi+,)Q—'/\Tii + l), 



(299) 



where, as usual, we introduced shifted arguments. If we denote H x with Q{xi,X2, xs, X4) the secondary 
quadrics and the corresponding coefficients are 

Q(rol) = Qi = Q(l -Xi,xi- X2, X2 - X3, X3 - X4), 

Q{T 12) = Q2 = Q{X1 -Xi,Xi- X2, X2 - X3,X3 - X4), 

Q{T2 3) = Q3 = Q{xi - Xi,X2 - X2,X2 - X3,X3 - X4), 

Q{T3 4) = Q4 = Q{xi - Xi,X2 - X2,X3 - X3,X3 - X4), 

Q{T = Q5 = Q{xi - Xi,X2 - X2,X3 - X3,-X^), (300) 

and Xq = 1,^5 = 0. Again, each coefficient in the Taylor expansion is written as a combination of 
threefold integrals which can be evaluated with standard BT techniques. This will introduce sub-leading 
quadrics, i.e. 

Qii = Qi{oi) = Qi{i,xi,x2), (301) 

etc, and sub-subleading quadrics, i.e. 



Qiji = Qij{0 1) = Qij{l,xi), 



(302) 



etc, and also constant terms, e.g. Qiji{l, 1) etc. At each step non-leading BT factors arc introduced and 
the procedure fails when one of the sub-leading BT factors is zero. In this case, since Eq is the sum of 5 
terms of the form 



1 f 

n=0 



with Qi constant, we rewrite the sum as 



,00 1 

^ n=0 ^ 



Likewise, by changing variables we obtain 



J dSsQiiQi + B^z)-^ = J dSsiX) dz qi {x^ Hi x + B^i + B5 z)'^ , 



(303) 



(304) 



(305) 



where B4i is the relevant sub-leading BT factor and we can apply a Mellin-Barnes splitting, followed by 
a sector decomposition, to x* HiX and B^i + Br,z^. For the Mellin-Barnes anti-transform the leading 
contribution comes from the pole at s = 3/2 giving 

(306) 



dz {B4i + B5 z)-^/2 ^ _ 
1- 



{B4^ + B,)-^' ~ B:'/' 



Alternatively we define Vo(xi, . . . , X4) = V{xi, . . . , 3:4) — B^ and write down the BT relations correspond- 
ing to 7^ and B^ = 0: 

1 



1 + (x - X) [Vo{xi, . . . ,X4) + B^]-^-'/^ = B5 [Vo{xi, . . . ,X4) + B^]-^-'/'', 



4 + e 



1 + 



1 

1T~e 



{x-x)d^ 



{xi, 



,X4) 



0. 



(307) 
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Then, after having summed each side of the two equations, we integrate by parts and set e = 0, obtaining: 

I r 4 

Eo = - dS2Y^ {X, - Xi+i) V-\i i + 1) Isub . (308) 

1=0 

Hence the study of the scalar pentagon reduces again to the study of four-point functions. The tensor 
pentagons too can be treated analogously to the case / 0. 

8.4 Evaluation of Eq when B5 ^ and Eq is singular 

For the pentagon we are in a special situation even if ~ and the condition Xi < Xi^i,i = 1, . . . , 5 
(with Xq = 1 and X5 = 0) is satisfied. Indeed, in this case, the BT algorithm already gives the correct 
leading behavior, Eq ~ I/-B5. This is easily seen by the Mellin-Barnes technique combined with a sector 
decomposition that would require to anti-transform the product 

r (.) r (3 - 3„, _ 1 
— 7:^2 ^5 ' p'-w 

Closing the s contour over the right-hand complex half-plane at infinity, thus selecting the poles at s = 2 
and s = 3 + k,k > 0, gives as the leading term. Therefore, in this case, we simply go on in treating 
the four-point functions according to their (sub-leading) BT factors. 



9 Infrared divergent Eq 

For the family, as long as -B5 7^ 0, we can use Eq. (|268| ) even for IR divergent configurations, therefore 
reducing the problem to the one of boxes. If instead -B5 ~ another procedure is more convenient. The 
physical relevance of the ii^-family in photonic virtual corrections to four-fermion production in e^e~ 
annihilation, as well as examples of the corresponding Eq reduction, have been illustrated in p^ j. 

Let us give an example for Eq where we consider the configuration with 

Pi=pf = — m^, mi = 0, m2 = = m, = m4 = M. (310) 

If we introduce the notation, 

[l]=g2, [2] = {q+pif +m\ [3] = {q + + p^f + M\ (311) 

[4] = (g + Pl+P2+P3)' + M^ [5] = (g-p5)2+m2, (312) 
then the following decomposition holds 

1 



jTirem 

Eq 



where the IR finite reminder is 



(pi +p2f + E'q^ + i {Dl^ + Dl% (313) 

-j^ y^^ [i]...[5] ' 

and where the two Dq functions are specified by their list of arguments: 

D^^ = Dq{pI ,pI,pI, {pi + p^f , {p2 + Pzf , (P3 + PA? ,m,M,M,m), 

Df = DQ{-m'^, {P2+P3f,pl, -m\ (P4+P5)^ {pi + p^f , , m , M , m). (315) 



47 



As done in Section the reminder is computed without simpHfications between numerator and denomi- 
nator. In particular, if 7^ we use the method of Section 8.1; otherwise, for = and Eq"^ regular 
(singular) we use the results of Section ^.3| (Section |8.4D . The classification of arbitrary IR divergent Eq 
functions follows from IR power counting similarly to the discussion of Section In all cases the IR 
divergent part is expressed through Cq functions that are contained in some of the Dq which appear in 
decompositions of the type expressed by Eq.(|313|). 



10 Six-point functions (F-family) 

To compute the generic, scalar, six-point function we introduce more compact notations: 



^^i- = (316) 

where Hpf is the (A^— 1) x (N—l) matrix with elements —pi-pj, = detHj^ is its determinant and ^N,ij 
is the co-determinant of the element Hij. The quadric will be written as = Hj^ x + 2 x + Lj^. 
It follows that 

b X 

Bn = = Gj^ — Kj^ Ajv A'jv, = '^N = ~ ■^jv Kf^. (317) 

Furthermore, = Tn\ and has elements 

KNi = \([i + l]-[i\), i = l,...,A-l, (318) 

where [1] = m\ and [j] = {pi -|- • • • -|- Pj-i)"^ + m'j, with j = 2, . . . ,N. Since the scalar hexagon is UV 
finite we can work in d = 4 dimensions where Gq = 0. Therefore, one BT iteration gives 



b6Fo = ^ I dSi Yl - ' (319) 

j=0 



where Xq = Xq = Q and where 



^60 


= y(o 1) 


= y(l,xi. 


• • ,2:4), 


^61 


= v{r2) 


= V{xi,xi 


• • • 


Vq2 


= V{2i) 


= V{xi,X2 


X2, ■■■ ,X4 


^65 


= V{fQ) 


= y(xi, •• 


,2:4,0), 



(320) 



and the scalar hexagon is the sum of six pentagons, i.e. of 30 boxes [24|. 

In computing the hexagon we have to decide about the input parameter set (hereafter IPS). Clearly 
internal masses, mi, • • • , mg belong to the IPS and for the rest we prefer to use invariants. For a total of 
N momenta we have 3 — 10 independent invariants which is considerably less than the number of inner 
products. The energy-momentum conservation leads to a number of linear relations among the inner 
products, of which a typical example is that connecting s, t, u for the box. Secondly, the dimensionality 
constraints lead to quadratic relations. For the hexagon there is just one, expressed by Gq = 0. The 
complete set of invariants will be (remember that in our conventions all momenta are flowing inwards) 
S'ijk- = — {Pi + Pj + Pk + • ■ O^, and we will select 8 of them for our IPS. 
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10.1 Form factors in the F-family 

Consider the integral 



1 /* Q ^ 

^'^ = i^^n^ I '^"'^ {q^+ml)--- {{q+pi+ ••• +^5)2 + ^2) = ^l^^^'^" ^^^^^ 
We immediately obtain 

j dS,V-'-^'\{x}) + Ft, (322) 



where F^- is a fourfold integral. Therefore there is a piece in proportional to -Fol*^ = 6), where the 
coefficient is 

5 

^ X,pi^. (323) 

't=i 

In order to prove that this coefficient is zero it is enough to contract it with a set of vectors pi that span 
d = 4 space-time. In this way we obtain 

5 5 
- Aij Kj pi-pi = Y^ Aij Kj Ha . (324) 

i=l i=l 

Since Ajj is the co-determinant of the element Hij we can use Laplace's theorem, Ajj Hu = 6ji G, to 
obtain 

5 

J2\jKjHa = GKi = 0, (325) 

1=1 

due to the vanishing of the Gram's determinant. For two or more powers of loop momenta in the 
numerator we need a new relation; let d = 4 — e and consider 

y 1 A °° 

x-;^)y'^-^ = — G = ^G(")e", (326) 

^ n=l 

where is arbitrary. We obtain 

XV''-^ = — -AdV, iG^^^fxV''-^ + G^^^XV^'-^ + —AdV^')= 0, (327) 
2 ^ 2 fj, 

etc. With the help of Eg. (|327|) we can show that the part proportional to Fq^cI = 8) in Fn^ has a coefficient 



5 



^ AijPi^Pj^, (328) 



which contracted with pi, I = 1, • • • , 5 gives 

5 



AijPi^Pj ■pi = - ^ Aij Hjipif, = - Gpif, = 0. (329) 

i,j=l i,j=l 
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In computing form factors of the F-family we will use results already obtained for the i?-family showing 
that, once more, everything is reducible to form factors of the L'-family. For instance for we use 



Xijpj^ = 0, — = —, Ge = 0, 

2=1 



and also 



4-e/2 ^ J_ 



1 + 



(x-X)d 



6 + e 

Furthermore, secondary quadrics are 



be 6 + e 



(330) 



(331) 



V{0 1) 


= V{l,Xi,X2,X3,X4) 


= ¥l(xi,X2,X3,X4), 


Vil2) 


= V{x, 


Xl 


X2 


X3 


xa, 


= V2{xi,X2,X^,Xa) 


V{2 3) 


= V{xi 


X2 


X2 


X3 


xa, 


= V3{xi,X2,X^,Xa) 


V{3l) 


= V{xi 


X2 


X3 


X3 


xa, 


= Va{xi,X2,X^,Xa) 


1/(45) 


= V{xi 


X2 


X3 


Xi 


xa, 


= V5(xi,a;2,a;3,a;4) 


V{f6) 


= V{xi 


X2 


X3 


X4 


0) 


= V(i{xi,X2,X3,XA), 



to which we apply a standard BT algorithm 



-3-e/2 



1 



1 + 



1 



4 + e 



or a special BT algorithm. 



1 



..-2-e/2 



(332) 



(333) 



(334) 



2 + e '""^ 

to show that is given by a combination of 813 threefold integrals, i.e. form factors of the Z)-family. 
10.2 Special configurations for Fq 



If we have = — A IT = on top of Gg = then Eg. ( 319 ) cannot be applied. In this case we can 
proceed as follows: since the vectors pi are linearly dependent, we substitute 



into the expression of V and show that 

5 4 



(335) 



HijXiXj - ^ H^j XiXj, Xi 

i,j=l i,j=l 



Xi 

Xi 77" X^^ 
^5 



(336) 



where H^^'i is a 4 x 4 matrix with elements — pi ■ pj . Likewise we have 

KiXi = KiXi + — X5. 
i=i 1=1 ^ 



(337) 
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If we define 

/ X \ 

Li = {l-ri)x5= \l- -^jx^, C/j = + (ri_i - ri)x5, (338) 
with xo = 1 and tq = 0, the expression for Fq can be cast into the following form: 

Fq = j dx^M^^), Mx5) = Yl I ^ dxi{Qi + 2^x^Y^ (339) 

where Q4 is a quadric in four variables, i.e. 

4 4 
Qa=Y^ H\f Xi Xj + 2Y^ Ki Xi + L. (340) 

i,j=l i=l 

To this quadric is associated a BT factor i?4 so that we can apply the BT algorithm to Eq.( |339| ) with a 
total BT factor + Ib^/ X5. 



11 Conclusions 

In this paper we presented a detailed investigation of the algorithms, based on the Bernstein-Tkachov 
theorem which form the basis for a fast and reliable numerical integration of one- loop multi-leg (up 
to six in this paper) diagrams. The rationale for this work is represented by the need of encompassing a 
number of problems that one encounters in assembling a calculation of some complicated process, e.g. full 
one- loop corrections to e^e~ 4fermions. Furthermore, in any attempt to compute physical observables 
at the two- loop level, we will have to include the one- loop part, and it is rather obvious that the two 
pieces should be treated on equal footing. 

Finally, our method represents a new strategy for the so-called problem of reduction of tensor integrals 
(for additional alternatives we refer to [p4|,[p^] and [p^). 

All algorithms that aim to compute Feynman diagrams numerically are based on some manipulation 
of the original integrands that brings the final answer into something smooth. This has the consequence 
of bringing the original (Landau) singularity of the diagram into some overall denominator and, usually, 
the method overestimates the singular behavior around some threshold. In these regions an alterna- 
tive derivation is needed. Instead of using the method of asymptotic expansions [^], we introduced a 
novel algorithm based on a Mellin-Barnes decomposition of the singular integrand, followed by a sector 
decomposition that allows us to write the Laurent expansion around threshold. 

Particular care has been devoted to analyze those situations where a sub-leading singularity may occur, 
and to properly account for those cases where the algorithm cannot be applied because the corresponding 
BT factor is zero although the singular point in parametric space does not belong to the integration 
domain. 

Clearly, no numerical evaluation of Feynman diagrams should be attempted unless we have a detailed 
knowledge of their analytic structure; in other words, we must know beforehand where the real singular- 
ities are sitting, which apparent singularity can be the origin of numerical instabilities, and what to do 
when one of these cases is met. 

One of the main by-products of the BT approach, Eq.(|l|), is that any diagram G has an associated 
factor Bq which is immediately computed for any one-loop diagram and gives information about the 
leading solution of the corresponding set of Landau equations (these solutions are notoriously hard to 
derive with standard methods |^^). Technically speaking, Bq = guarantees that the Landau equations 
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admit a proper solution. Repeated applications of the algorithm, as explained for instance in Section jA. , 
will also introduce sub-leading Sgig^-factors whose zeros correspond to sub-leading Landau singularities. 
Likewise, our solution around Bq = will introduce sub-leading ^-coefficients of a second kind (as 
explained, for instance, in the discussion after Eq.(^)), not associated with Landau singularities of the 
reduced diagrams. 

When these pieces of knowledge are at hand, we can automatize our calculations, even without 
having to worry in advance about the physical nature |31] of the singular points (a much more difficult 
assignment): it will be enough to be able to handle all cases and let the program decide when Bq = 
corresponds to a singular behavior or when a coincidence occurs. Furthermore, we have a conjecture 
stating that when Bq = 0, then the simultaneous occurrence of Bq^^^ = is not physical, but we have 
not been able to prove it. Instead we decided to cover also this case, so that our expressions for the 
one- loop diagrams are quite general and not limited to the physical region of a specific process. 

In general it is known that when a given Landau singularity curve touches one of its sub-leading 
singularity curves the determination of the nature of the singularity breaks down completely (because 
the ai's are no longer uniquely determined) and changes suddenly and drastically. Similarly, our analysis 
via the multiple Mellin - Barnes techniques (cf. Eqs.(p^)-(p3|) and Eq.( |221| )), dealing with the coincidence 
Ba = = (where B^^g^ is of second kind), represents a novel result. Alternatively, we have been 

able to derive new integral representations that bypass the use of Laurent expansions. 

Finally we have given a description of infrared divergent one-loop virtual configurations in the frame- 
work of dimensional regularization: here both the residue of the infrared pole and the infrared finite 
remainder are cast into a form that can be safely computed numerically. 

The collection of formulas that cover all corners of phase space have been translated into a set of 
FORM codes [^] and the output has been used to create a FORTRAN code whose technical description 
will be given elsewhere. 
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A A master integral 

In this appendix we calculate the master integral 



Jo 



(341) 



where a and h are real numbers with h approaching and a far from 0. If /? > 2a — 1 the integral is 
not divergent, otherwise we split the integral into two pieces: the first one contains the divergent part as 
6 — >• 0, while the second is finite: 



I{a,i3)ia,b) = J dxx'^ {ax"^ + b - i5)~" - J dx x^ {ax^ + b - id)'"" . 



(342) 



Since we are in the limit of 5 ^ 0, we may write ax"^ + b — i5 = {a — i5)x'^ + b — id. Next we perform the 
following transformations in the two integrals respectively: x — > x^^^ and x —> x~^ and obtain 



<3-l 



/(cf,/3)(a, &) = - / dxx 2 (^a — iS) X + b — i6 



dx X 



a — i6 + {b — i6) x^ 



In the first integral an additional change of variable is performed: 

b — id 1 — X 



a — i6 X ' 



and we make use of the following property: 



b-i5Y {b-i5Y 



a — i5 J {a — i5)^ 



Hence the final result is: 



/3+1 



^(a,/3)(a,&) 



2{a-i5)^ 



B 



-, a 



a — i5 + {b — i5) x^ 



(343) 



(344) 



(345) 



(346) 



The second term is well defined only for /3 < 2q — 1 and we perform an analytic continuation to the 
region (3 <2a + \ (with a,b ^ a,b — i6): 



I{a,f3)ia,b) 



2a 2 



B 



, a 



2a- (3-1 



dx X 



2a~l3 



\a + bx 



(347) 



From Eq.(347) we see that the two divergent terms cancel out also when (3 = 2a — 1. For a = n and 
/3 = 2 n — 1 we get 



'(",/3) 



(a, 6) 



dx 



[a + bx^] 



2l-n 



(348) 



53 




54 



References 



G. 't Hooft and M. Veltman, Nucl. Phys. B 153 (1979) 365. 

A.V. Kotikov, Phys. Lett. B 254 (1991) 158; 
A.V. Kotikov, Phys. Lett. B 267 (1991) 123. 

L.D. Landau, Nucl. Phys. 13 (1959) 181. 

R.J. Eden, P.V. Landshoff, D.I. Olive, and J.C. Polkinghorne, The Analytic S-Matrix, Cambridge 
Univ. Press, 1966. 

W.L. van Neerven and J. A. Vermaseren, Phys. Lett. B 137 (1984) 241; 

Z. Bern, L.J. Dixon and D.A. Kosower, Nucl. Phys. B 412 (1994) 751 JarXiv:hep-ph/9306240| ; 
A.I. Davydychev, J. Math. Phys. 33 (1992) 358. 

A. Vicini, talk given at Loopfest, Brookhaven National Laboratory, Upton, NY, May 9-10, 2002 
( http: / /www.pas.rochester .edu^dow /loopfest / parallel.html ) . 

W. Beenakker et al, Nucl. Phys. B 500 (1997) 255 ||arXiv:hep-ph/961226q] . 

G. Passarino and M. Veltman, Nucl. Phys. B 160 (1979) 151; 

D. Bardin and G. Passarino, The standard model in the making: Precision study of the electroweak 
interactions, Oxford, UK: Clarendon (1999). 

G. Passarino, Nucl. Phys. B 619 (2001) 257 [|arXiv:hep-ph/0108252||; 

G. Passarino and S. Uccirati, Nucl. Phys. B 629 (2002) 97 |arXiv:hep-ph/0112004 |. 



D.Y. Bardin, L.V. Kalinovskaya and F.V. Tkachov, |hep-ph/0012209 



G. Passarino, |hep-ph/0101299 



F.V. Tkachov, Nucl. lustrum. Meth. A 389 (1997) 309 ||hep-ph/960942g| ]; 
L.N. Bertstein, Functional Analysis and its Applications, 6(1972)66. 

A. Erdelyi et al.. Higher Transcendental Functions, vol. 2, Bateman Manuscript Project, McGraw- 
Hill, 1953; 

L.J. Slater, Generalized Hypergeometric Functions, Cambridge Univ. Press, 1966. 

T. Binoth and G. Heinrich, Nucl. Phys. B 585 (2000) 7 41 [|arXiv:hep-ph/00040"Tl ; 
M. Roth and A. Denner, Nucl. Phys. B 479 (1996) 495 [|arXiv:hep-ph/9605420|| ; 
K. Hepp, Commun. Math. Phys. 2 (1966) 301. 

J. Fujimoto, Y. Shimizu, K. Kato and Y. Oyanagi, Prog. Theor. Phys. 87 (1992) 1233; 
J. Fujimoto, Y. Shimizu, K. Kato, N. Nakazawa and T. Kaneko, prepared for 10th International 
Workshop on High-energy Physics and Quantum Field Theory (NPI MSU 95), Zvenigorod, Russia, 
19-26 Sep 1995; 

F. Caravaglios, Nucl. Phys. B 589 (2000) 475 [|arXiv:hep-ph/0004030|1 ; 



R. Easther, G. Guralnik and S. Hahn, Phy s. Rev. D 61 (2000) 1 25001 ||arXiv:hep-ph/9903255|| ; 
M. Steinhauser, Phys. Rept. 364 (2002) 247 1arXiv:hep-ph/0201075i . 



55 



[15] G. Montagna, F. Piccinini, O. Nicrosini, G. Passarino and R. Pittau, Comput. Phys. Commun. 76 
(1993) 328; 

G. Montagna, F. Piccinini, O. Nicrosini, G. Passarino and R. Pittau, Nucl. Phys. B 401 (1993) 3. 

[16] E.E. Boos and A.I. Davydychev, Theor. Math. Phys. 89 (1991) 1052 [Teor. Mat. Fiz. 89 (1991) 56]; 
A.I. Davydychev and J.B. Tausk, Nucl. Phys. B 397 (1993) 123; 



S. Bauberger, F.A. Berends, M. Bohm and M. Buza, Nucl. Phys. B 434 (1995) 383 |arXiv:hep- 
ph/9409388i; 



A.I. Davydychev and J.B. Tausk, Phys. Rev. D 53 (1996) 7381 |arXiv:hep-ph/9504431 |. 

[17] R. Karplus, C. M. Sommerfield and E. H. Wichmann, Phys. Rev. Ill (1958) 1187; 

G. Kallen and A. S. Wightman, K. Dan. Vidensk. Selsk. Mat. - Fys. Skr. 1 NO. 6 (1958). 

[18] R. E. Cutkosky, J. Math. Phys. 1 (1960) 429. 

[19] G. Devaraj and R.G. Stuart, Nucl. Phys. B 519 (1998) 483 |arXiv:hep-ph/9704308(l . 

[20] L. Lewin, Poly-logarithms and Associated Functions, North Holland (New York 1981). 

[21] W. Beenakker and A. Denner, Nucl. Phys. B 338 (1990) 349. 

[22] H.D. Abarbanel and C. Itzykson, Phys. Rev. Lett. 23 (1969) 53; 
M. Levy and J. Sucher, Phys. Rev. 186 (1969) 1656. 

[23] G. Sterman, An Introduction to Quantum Field Theory, Cambridge Univ. Press, 1993. 



[24] J. Fleischer, F. Jegerlehner and O.V. Tarasov, Nucl. Phys. B 566 (2000) 423 iarXiv:hep-ph/9907327i ; 
O. V. Tarasov, Acta Phys. Polon. B 29 (1998) 2655 ||arXiv:hep-ph/981225q] ; 
Z. Bern, L.J. Dixon and D.A. Kosower, Nucl. Phys. B 412 (1994) 751 ||arXiv:hep-ph/9306240 l. 



[25] T. Binoth, J.P. Guillet and G. Heinrich, Nucl. Phys. B 572 (2000) 361 [|arXiv:hep-ph/9911342| ; 
Z. Bern, L. J. Dixon and D. A. Kosower, Phys. Rev. Lett. 70 (1993) 2677 [ |arXiv:hep-ph/930228q| . 

[26] A. Denner, S. Dittmaier and M. Roth, Nucl. Phys. B 519 (1998) 39 iarXiv:hep-ph/9710521 1. 

[27] T. Binoth, J.P. Guillet, G. Heinrich and C. Schubert, Nucl. Phys. B 615 (2001) 385 ||arXiv:hep- 



ph/0106243 |. 



[28] C. Bauer and H.S. Do, Comput. Phys. Commun. 144 (2002) 154 ||arXiv:hep-ph/0102231 



S. Weinzierl, Phys. Lett. B 450 (1999) 234 iarXiv:hep-ph/9811365| ~ 
A.I. Davydychev, Phys. Lett. B 263 (1991) 107. 

[29] F.V. Tkachov, in Quarks-82 Proc. Int. Seminar, Eds: A.N. Tavkhehdze et al. INR, the USSR 
Acad. Sci. Moscow, 1982; Euclidean asymptotic of Feynman integrals, preprint INR P-332 (Moscow 
1984); Int. J. Mod. Phys. A8 (1993) 2047; Sov. J. Part. Nuclei 25 (1994) 649. 

[30] A. C. Wu, Phys. Rev. D 9 (1974) 370. 

[31] D. Kershaw, Phys. Rev. D 5 (1972) 1976. 

[32] J.A. Vermaseren, |niath-ph/0010025| . 



56 



